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Abstract 

The most general basis of operators parametrizing a low-scale departure from the SM particle 
content is constructed. The SM gauge invariance is enforced, and operators of lowest dimensions 
are retained separately for a new light neutral particle of spin 0, 1/2, 1, and 3/2. The basis is 
further decomposed into couplings to the SM Higgs/gauge fields, to pairs of quark/lepton fields, 
and to baryon/lepton number violating combinations of fermion fields. This basis is then used to 
systematically investigate the discovery potential of the rare FCNC decays of the K and B mesons 
with missing energy in the final state. The most sensitive decay modes in the s — > d, 6 — > d, and 
6 — > s sectors are identified and compared for each type of couplings to the new invisible state. 
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1 Introduction 



Though extremely successful, the Standard Model (SM) is not expected to be valid up to arbitrary high 
energies. It certainly needs to be amended at the Planck scale, with the advent of quantum gravity. But 
then, some new dynamics should show up already at a much lower scale, though still above the electroweak 
scale, to avoid hierarchy problems. With this picture in mind, it seems natural to assume that all the New 
Physics (NP) degrees of freedom are heavier than the known SM particles, and decouple at low energy. In 
a fully model-independent way, the impact on the SM of any NP model can then be embedded into non- 
renormalizable effective operators, under the provision that these involve only the SM fields and satisfy the SM 
gauge symmetries. The full classification of these operators has been achieved many years ago, by Buchmuller 
and Wyler [I] for the baryon- and lepton-number conserving operators, and by Weinberg [2] for those violating 
these global symmetries. 

Still, whether the SM particles are the only dynamical degrees of freedom within the electroweak energy 
range is far from established. Not only is the existence of new light particles not excluded, since they would 
evade direct detection when sufficiently weakly interacting, but their presence could even be welcome. Indeed, 
many NP models are built upon some spontaneously broken symmetries, and do often have remnants at low- 
energy in the form of massless or very light Goldstone bosons. A well-known example is the axion [3 , introduced 
to cure the strong CP problem of the SM. More crucially, there are now very strong indications that the universe 
is filled with dark matter, so there should be at least one new electrically neutral colorless particle, possibly 
lighter than the electroweak scale. Once opening that door, it is not such a drastic step to imagine a whole 
dark sector, i.e. a full-fledged set of darkly interacting dark particles only loosely connected to our own visible 
sector. Further, it should be stressed that adjoining a dark sector to the SM is always possible, does not need to 
be directly related to dark matter (so one would rather speak of a hidden sector) , and is actually quite generic 
in supersymmetric models. For a recent review, including further physical motivations from string theory or 
extra dimensional settings, see e.g. Ref. H]. 

In the present work, our main goal is to construct the lowest-dimensional effective interactions parametrizing 
a low-scale departure from the SM particle content. Specifically, we assume that there is a new particle of spin 
0, 1/2, 1, or 3/2, neutral under the SM gauge group SU{3)c SU{2)l (Xi U{1)y, and write down the gauge 
invariant operators coupling this particle to SM fields. These effective interactions are not yet included and 
thus complement the NP operator basis of Refs. [Tl[2]. To our knowledge, such a complete basis has never been 
presented, though parts of it already appeared in the literature. In particular, those effective couplings between 
SM and dark particles which are renormalizable, sometimes called portals, have already been investigated [S]. 

Our second goal is to constrain the effective operators. Since the new state is assumed neutral under the 
SM gauge group, it looks like a natural dark matter candidate. However, the viability of this hypothesis 
would require constraining its mass, couplings, and lifetime, and this in general requires more inputs about 
the dark sector dynamics. Here, we refrain from doing so and rather concentrate exclusively on the quark 
FCNC transitions s — > dX, b — >■ dX, and b — > sX, where X collectively denotes any final state made of dark 
particles. Our focus on these modes, instead of for example leptonic observables or collider signals, is motivated 
on one hand by their extreme sensitivities to NP (as detailed in the next section), and on the other, by the 
next generation of experiments currently under construction. Indeed, rare K decays are the main targets of 
the NA62 (CERN) and KOTO (J-Parc) experiments, while rare B decays could be accessed at the Super-B 
(Italy) and Belle II (KEK) facilities. So, in the present work, we further assume that the dark particle is light 
enough to be directly produced in K and/or B meson decays, and sufhciently long-lived to escape detection in 
flavor factories. For all practical matters, this new particle is invisible, and would show up as missing energy 
in FCNC-induced rare K and B decays (for recent works along this line, see e.g. Refs. [BIH)- 

Our analysis is organized according to the dark particle spin. To keep the discussion focused on the operator 
basis, we rely on extensive appendices to cover the issues of hadronic matrix elements and differential rates. 
So, before entering the discussion, the next section summarizes the main features of the observables considered 
here, i.e. the rare K and B decays. Then, given our focus on FCNC processes, a flavor-based classification of 
the dark operators is described in the following section, on which we rely throughout the paper. 
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Figure 1: The rare decays with missing energy in the SM, as induced by the Z penguin (W boxes are understood). 
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Table 1: Naive reach, in terms of scales A and as a function of the effective operator dimension n, of the rare 
FCNC-induced K and B decays, as estimated from Eq. ^ with the CKM values of Eq. ©. 

Rare FCNC decays 

The FCNC-induced decay modes are very suppressed in the SM, where the missing energy is carried away by 
a I'D pair (see Fig. [1]). So, even relatively small NP contributions could be evidenced. Specifically, to set the 
stage and get an idea of the sensitivity of the rare K and B decays, imagine that a NP operator of dimension 
n contributes to d^ — >■ d"^ X, with 7 = 2,3, J = 1, 2 the quark generation indices. If its Wilson coeflacient is set 
to one, then there is a scale A such that the NP contribution equates the SM prediction for d^ — >■ d'^ vv, 

A^^ m1;16^'^**^^''^'' 

with TO2,3 — ^K,B, 9 the SU{2)l coupling constant, and V the CKM matrix. As shown in Table [1] as a 
function of the dimension, the SM loop factor combined with the CKM suppression pushes the scales A well 
above the electroweak scale for n < 7. On the contrary, the rare decay constraints cease to make sense for 
SU(2)l <8) U{1)y invariant operators of dimension n > 9, since powers of [H'^ H)/ h? — >■ / h? grow unchecked 
when A < where v « 246 GeV is the SM Higgs vacuum expectation value. 

Clearly, these scales are only indicative. The true sensitivity to a given dark operator depends essentially on 
two additional factors. First, the quark transitions d^ d'^ X have to be probed through hadronic processes. 
Hence, depending on the modes, hadronic matrix elements as well as phase-space factors can alter significantly 
the estimates of Table [T] In the following, we compare the sensitivities of all the leading modes. Specifically, in 
the K sector, we include the modes with the least number of pions and photons in the final states, i.e. K X^ 
K ttX, K -fX, and K ttttX, and leave out the K — ttttttX modes. Similarly, in the B sector, the 
considered modes are the B ^ X , B ^ {K, K*)X, and B — > (tt, p)X decay channels. The jX channel, driven 
in the K sector by the QED anomaly, is suppressed and difficult to reconstruct experimentally in the B sector, 
and will thus not be included [B]. 

The second factor determining the true sensitivity of a given mode is related to the experimental strategies 
deployed to measure it. Since invisible states are not seen, the kinematical reconstruction is limited. In addition, 
these modes are so rare (in the SM) that they require very aggressive background suppressions. To this end, 
the central tool is the differential rate in terms of the kinematical parameters of the visible products. But this 
differential rate depends on the nature of the dark particle. Currently, most experimental analyses implicitly 
impose the SM differential rate (for X — uD) . This means that the current bounds cannot be directly translated 
to other types of final state particles. This motivates another goal of the present paper, which is to provide the 
full dictionary of the differential rates for all the leading effective interactions involving invisible final states. 
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These spectra should be used by the experimentahsts to derive bounds for each type of new invisible state. 

For more details on these issues, including kinematics, matrix elements, current measurements or bounds, 
and experimental prospects for the various modes, we refer to appendix I A. II (IB . 1 p for K (B) decays. 

Flavor-based classification of the dark operators 

At the electroweak scale, once the whole NP particle spectrum but the X has been integrated out, the lowest 
dimensional operators can be split into three types according to their quark and Icpton field contents: 

Tieff = 'Hmat + 'Hint + 'H^B-AC ■ (2) 

By definition, "Hint contains only gauge and Higgs fields, while Umat and Habac contain at least one SM 
fermionic field. The operators of Hab.ac have a non-zero charge under the baryon (B) or lepton (£) num- 
ber U{l)s. As in Ref. 'Ij 2i, all the operators are to be written as manifestly invariant under SU{3)c <E) 
SU{2)l (X" U{1)y, i-C. in terms of the quark (lepton) doublets Q (L) and singlets U,D (E) of each flavor, 
of the SU{3)c, SU{2)l,U{1)y field strengths G^^,, PV^^, 5^,^, of the Higgs doublet H, as well as of covariant 
derivatives acting on these fields, insofar as these cannot be reduced using the SM equations of motion (EOM). 

Due to their different field contents, these three types of operators do not contribute equally to the quark 
transitions s — > dX, b — >■ dX, and b — > sX, so let us organize them differently, in terms of four classes of 
scenarios, as shown in Tabled 

Consider first the operators of Umat involving down-type quarks (those with leptons or up-type quarks are 
obtained by substituting D,Q E, L oi D, H U,H*). Up to possible partial derivatives acting on the 
quark or invisible fields, and omitting the Dirac structures, the quark currents are 

-Hmat = ^hWq' xX+ ^HQ'D^ xX + ^Q'Q^ xX+ ^d'D-' xX . (3) 
A" A" A" A" 

Those operators have a flavor structure, and thus can in principle induce d^ — > d'' X. Clearly, when analyzing 
the physics reach of rare K and B decays in terms of the scale A, the assumptions made on the c^^'^ are crucial. 
There are two main scenarios: 

I. The constraints derived from rare FCNC decays are the tightest when the NP flavor structure is generic, 

c'^J ^ 0(1) . (4) 

As shown in Table [l] the bounds on the NP scale A are then often far above the electroweak scale. 

II. Since T-Lmat results from integrating out the whole NP particle spectrum, the flavor-breaking character 
of its operators could originate from dynamical effects not related to the dark sector. In that case, the 
NP dynamics would also quite naturally correct the visible FCNC operators, on which there are many 
tight experimental constraints from K and B physics [5]. This is typically the case in supersymmetric 
settings, where the flavored soft-breaking terms cannot be fully generic. Phenomenologically, a simple 
way to account for such a non-generic NP flavor structure is to impose the Minimal Flavor Violation 
(MFV) ansatz [5], i.e. force the quark currents to have the forms 

D\YaYlYuY'Q' , Q^(YtY„)"Q^ , D\Y ^YIY ^Y^^' ' D' . (5) 

In the down-quark mass-eigenstate basis, the diagonal vYd — ^/2va.d tunes the chirality flips, while 
vYu = \/2u\uV parametrizes the flavor change (m„_rf denotes the diagonal quark mass matrices, V the 
CKM matrix, and v 246 GeV the Higgs vacuum expectation value). So, MFV rescales the Wilson 
coefficients according to c^;^ = m^c£'^/t), c£';[j — c^j^j^m'^v, c^j^j^ — m^m^c^'^/z;^, and 

r X'"^ (-3.1 + ^1.3) X 10-" , 
- = Yt Y„ « v;,Vtj X"'- (7.8 - i3.1) X IQ-^ , (6) 

[ A''* w (-4.0 - i0.07) X 10-2 . 

Upon these rescalings, the accessible scales A are then much lower, especially for operators of low dimen- 
sions, and for s — >■ d operators involving light right-handed quarks. 
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If the whole NP dynamics is flavor blind, then X couples only to the flavor-diagonal quark currents and 
c^^'' = 0. A flavor transition is of course still possible but it must proceed through the SM weak interactions, 
i.e. the flavor-blind quark currents must be dressed by a flavor-changing W interaction, either between the 
quark lines or as a self-energy on these quark lines. 

In particular, all the operators of %int are of this type. Indeed, to contribute to FCNC processes, gauge 
fields have to be coupled to quarks, while Higgs fields are either coupled to quarks or left as external tadpoles 
(to be replaced by the vacuum expectation value after the electroweak symmetry breaking). The resulting 
couplings between quarks and X can then be matched onto Hmat, and satisfjQ c^^'-' « 0. Thus, the only 
difference with respect to the flavor-blind Hmat operators is that the c^^ coefficients are initially suppressed 
by some power of the SM coupling constants (at the scale A), by some loop factors, and by quark Yukawa 
couplings when a Higgs field is coupled to the quark line (so that c^^ ^ c^^'^^). 

For reasons entirely pertaining to the SM dynamics, it is different to probe the couplings to heavy quarks, 
c^^, and to light quarks, c^^'^^, so these constitute our third and fourth classes of scenarios. 

III. Let us assume we have an operator coupling X to the top quark current. Dressing it with a W exchange 
(see Table[2]), the necessary GIM breaking arises at the electroweak scale. From the rare decay perspective, 
this electroweak physics is local, so it can be matched onto the flavor-changing operators of Class II. The 
Wilson coefficients end up suppressed by the MFV scalings ([5]) and by a loop and gauge coupling, i.e. 



Typically, the bounds on the scale A are brought down very signiflcantly, often at around the electroweak 
scale. Still, the rare K and B decays remain ideal probes for such kind of effective couplings since a direct 
collider signal in the ti channel is presumably hidden by the large flavor-blind SM backgrounds. 

IV. With X coupled to the light quarks (u, d, or s), it is much trickier to derive bounds on the scale A for three 
reasons. Firstly, B physics is unable to provide competitive constraints, given the small CKM factors 
V*i,Vus or V*f^Vud- Secondly, though in the K sector the weak transition is favored by the large V*gVud 
for CP-conserving observablef@, the light quarks are never integrated out and remain dynamical. This 
renders the theoretical control difficult since the K physics scale is too low to allow for a perturbative 
QCD treatment. This is true both for the effective operators describing the weak transition and for 
the effective operators describing the production of the new invisible states. Thirdly, there are already 
many constraints on the flavor blind production of invisible particles, in particular from 7r° or quarkonium 
decays. Compared to the other classes, it is a priori not clear whether rare decays offer privileged windows. 
In the following, we will consider only speciflc scenarios where competitive bounds can be derived. 

The flnal type of operators is a bit different and does not immediately flt in the above classification. First, 
in general, AB and AC dark operators cannot be present simultaneously, since an X exchange would induce 
proton decay. A way out would be to impose MFV, which forces AC operators to be proportional to the 
tiny neutrino masses [IT], rendering them irrelevant. The AB operators are not particularly suppressed under 
MFV, but they never contribute to the FCNC-induced rare decays considered here. Indeed, these modes 
trivially conserve B since an odd number of baryons would be required in the final state. This is not possible 
for K decays, while the signatures for B decays should be experimentally clear, but are beyond our scope. 

If MFV is not imposed and AB operators are somehow disposed of, then AC contributions to the — >■ d"^ X 
processes are possible. But, the only AC ^ invisible final states either involve an odd number of neutrinos, 

^The contributions to c^^^ are tiny because, as long as the SU{2)l ® U(1)y symmetry is exact, the gauge/Higgs fields of the 
Hint operators can couple to a flavor-changing quark current only at the cost of at least two Yukawa insertions, i.e. two Higgs 
tadpoles, on the quark line. At the high scale, these tadpoles cost a factor A"'^ compared to the flavor-blind version of the same 
operator. When A is large, flavor-changing effects thus dominantly originate from electroweak scale GIM breaking, even though 
this necessitates an extra W loop (see Table [2jl. 

^CP-violating observables, for which Im(yj'gV„d) = (this standard CKM convention is used throughout the paper, see e.g. 
Ref. |10p. are induced by heavy quarks and fall into Class III. Due to the scaling (O, they are very suppressed. 




(7) 
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Table 2: Flavor-based classification of the operators involving dark particles, collectively denoted X. After the 
electroweak symmetry breaking, the Hmat operators are directly matched onto 7ieff{q^ — > Q'^^) and split into 
the four classes. The Hint operators collapse onto the Class III and/or IV flavor-blind operators once their 
gauge/Higgs fields are attached to quarks. The Hab,ac operators have different signatures, and do not fit in 
this classification. Note that the charm quark is considered heavy (light) for K (B) decays. 

or some AC — 2 neutrino pairs vl^l- Since a neutrino field in an effective operator costs A^'^/^, these are in 
general significantly suppressed compared to the operators of Hmat and Hint ■ The only exceptions are those 
contributing to — )■ d^VL4' or d^ — )■ d^VL'^. As will be discussed in the relevant sections, because vl is part 
of the lepton doublet, these operators are always accompanied by the charge-current transitions d^ — > u'^i~ip 
or d^ — > u'^i^'i', which may offer better windows than the rare FCNC transitions. 



2 Invisible spin-1/2 fermion 



When the new invisible fermion is neutral under the SM gauge group and is produced in pairs, imposing the 
SU{3)c fX" SU{2)l ® U{1)y gauge invariance reduces the basis to the usual ten chiral currents: 

s s s s 

+ ^H^DQ X i^Li'R + ^HWQ X i,RtPL + ^HQD x V^lV'a + ^HQD x i/i^VL 

T T 

+ '^H^Da^M X V^fla^'^^L + ^HQa^^D x ^icx'^^i^ , (8) 

with the definition a^^^ = i{'-f'^^'^ + g'"^), and where Q (D) stands for the left-handed quark doublet (right- 
handed down-quark singlet). Similar operators can be written down for the up-quark right-handed singlet or 
for leptonic transitions, and the generalization to a two Higgs-doublet model is straightforward. 

The coefficients are not assumed real, and their flavor indices are understood. For example, written in full 
for the s — >■ d, 6 s, and b ^ d sectors which concern us here: 

V V,sd V,bs V,bd 

'^Qlt.Qxi^Ll^i^L = ^^s-t^,PLdxi)Ll^^i)L + '^h^.PLSx4,Ll^^PL + ^j^h^.PLdxi)Ll^iJL + h.c. . (9) 
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Table 3: Pair production of invisible states: Scales A (in TeV) accessible for each operator, assuming 10^ 
bounds on all the branching ratios. The operator dimensions d are written with a bar for those involving a 
Higgs field, and thus aw/A factor. For the K — )■ ttXX modes, the differential rates are integrated over the 
pion momentum window [140, 195] U [211, 229] MeV, see App. IXTl The short-hands XS (XP) stand for SS or 
PS (SP or PP), as appropriate. The quoted m'^'^^ are indicative; the bounds are derived for mx = or, if 
the rate vanishes (signaled by (*)), for mx = 100 MeV. For those which do not vanish, most V — A, S ~ P, 
01 T — T degeneracies are lifted when mx 7^ 0. Still, the dependences of the As on mx are rather weak. The 
bounds stay within an order of magnitude of their values at mx = 0, except for mx close to the kinematical 
threshold, where they sharply drop towards zero. This is shown in the last column for — > tt+i/i/, where A 
(normalized to the value quoted in the table) is plotted against mx over the range [0,rn7r]- 



There are only two tensor operators because the identity 2a^'^^z — ie'^^'^"aap forbids Qa^^D x ^pfja'^'^ipL and 
DcTiiuQ X ipLcr'^'^ipR- The dimension-seven operators involving the SM Higgs field are retained because their 
suppression is only by v/ A. By contrast, true dimension-seven operators involving an additional derivative are 
more severely suppressed by mK.B,i>/A, and are thus not included. 

After the electroweak symmetry breaking, the gauge-invariant operators are rewritten in terms of the 
(pseudo)scalar, (axial)vector, and (pseudo)tensor currents. This minimizes the interferences between the cur- 
rents in physical observables, since these operators dominantly produce the invisible fermions in different states. 



7 





d 


^ XX 


B nXX 


B pXX 


B KXX 


B K*XX 


^max 




mgo /2 


{■niB - m^)/2 


{niB - mp)/2 


{niB - mK)/2 


{ruB - ■mK*)/2 


1/2 


fss.sp 


7 


— 


3.2 


— 


3.0(4.6) 


— 




fps,pp 


7 


2.7 


— 


2.1 


— 


1.7(4.0) 




fvV,VA 


6 


- 


10 


5.3 


9.7(19) 


4.3(15) 




Iav 


6 


- 


- 


7.6 


— 


6.2(21) 




fAA 


6 


6.2* 


— 


7.6 


— 


6.2(21) 




Itt 


7 


— 


2.6 


3.1 


2.4(3.7) 


2.7(6.2) 




ffr 


7 


— 


2.6 


3.0 


2.4(3.7) 


2.6(6.0) 





gss 


6 


— 


82 


— 


75(140) 


— 




gps 


6 


50 


— 


46 


— 


37(130) 




9VV 


6 


— 


7.2 


— 


6.8(13) 


— 




9AV 


6 


— 


— 


4.7 


— 


3.9(13) 


1 


hss,sp 


8 


— 


0.61 


— 


0.56(0.79) 


— 




hps,pp 


8 


0.56 


— 


0.42 


— 


0.36(0.67) 


3/2 


fss,sp 


9 


— 


0.18 


— 


0.17(0.22) 


— 




fps.pp 


9 


0.18 


— 


0.13 


— 


0.10(0.19) 




fvV,VA 


8 




0.15 


0.10 


0.14(0.19) 


0.09(0.17) 




fAV 


8 






0.12 




0.11(0.20) 




fAA 


8 


0.14* 




0.12 




0.11(0.20) 






9 




0.16 


0.16 


0.15(0.19) 


0.15(0.24) 




fTS,TP 


9 




0.14 


0.12 


0.13(0.17) 


0.11(0.18) 




ffs 


9 




0.11* 


0.14 


0.09* (0.12*) 


0.12(0.20) 




ffp 


9 




0.13* 


0.14 


0.11* (0.15*) 


0.12(0.20) 



Table 4: Pair production of invisible states: Scales A (in TeV) accessible for the various operators with present 
(future) measurements (see App. [BTTI) . The operator dimensions d are written with a bar for those involving a 
Higgs field, and thus a. v/A factor. We assume mx = everywhere, except when the rate vanishes in this limit 
(indicated by *). Specifically, the B — >• XX bound on fAA, the B — ^ ■jt{K)XX bounds on gr, fTS,TP, and the 
B — >■ K*XX bounds on gT,v,A are derived for mx — 2 GeV. The behaviors of the scales A as functions of mx 
are shown for X = in Fig. [51 
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and flavor indices are under- 



where fxY tunes q^Txq'^ x iprY4' with F 
stood. The corresponding differential rates are listed in Appendix lA. 31 (jB.3l) for K (B) decays, and the physics 
reach is summarized in Tables [3] . Specifically, the entries of these tables are obtained by turning on each 
fi in turn, while keeping the others to zero. We do not do this at the level of the Ci in order to minimize the 
interferences among the NP contributions. The value of each |/,;| is set to one for the vector and axial vector 
currents, and to v/A for the others, in order to keep track of the SU{2)l €5 U{1)y structure of the underlying 
operators. Since K^^s are approximate CP eigenstates, the CP phase of the fi must be kept arbitrary. Effec- 
tively, we turn on Im fi and Re fi to one (or v/A) separately, and the tightest bound (largest A) is indicated 
in Table H 

To derive the entries of Tables [3] and SI specific experimental bounds on the rare decay branching ratios are 

for K (B) decays. Note that the scales A corresponding to tighter or looser 



used, as detailed in App. lA.ll (IB.ip 
experimental bounds are immediately obtained by a simple rescaling of the numbers in Tables |3] and |4j given 
the dimensions of the operators ([8]) and the definitions (ITU)) . 
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Figure 2: Evolution of the scales A associated with the pair production of invisible fermions as a function of 
the ferniion mass. The values at = correspond to those quoted in Table El except for B — )■ ip-ip, given 
at = 2 GeV. Note that these B {K^*\ p^Tr)'ipip bounds assume flat experimental acceptances and full 
phase-space coverage, which is not true in the present experimental analyses, see App. IB.ll 



Finally, in Fig. |31 we compare the sensitivity of the various K and B decay modes for two illustrative 
examples, c^^ and cf^, with and without imposing MFV. In the former case, the coefficients are set to 
{cY^lY'^ = A^"^ and (cf^)^"^ — \^''mdi/v, see Eq. Actually, to avoid dragging a relative factor nib/v, the 
scales A are plotted taking a two Higgs doublet model of type II at large tan/3 = Vu/vd (where u„ and Vd are 
the two Higgs vacuum expectation values), so that nib/vd ~ 1 and (cf^)'*'^ is simply suppressed by nis/mi,. As 
can be seen in Fig. [31 this chirality flip is expensive in the K sector, and pulls the K TTipip constraint an 
order of magnitude below those from rare B decays. By contrast, for the vector current, similar constraints 
are drawn from rare K and B decays when MFV is imposed. 

Let us now turn to the operators not involving the SM fermions, but rather the SM gauge and Higgs fields. 
Using these fields EOM, as well as partial integration and identities like 2a'^^j5 = ie'^'^^^aap, the leading 
operators reduce to 



+ ^iH^^^H X 7/ii7^Vi + ^iH^^^H X ^R-ff'^R 

pG pG 
, ^RL,i 7 , , ^LR.,i j , 

'-Qi X VflV'L + — pj^Qi X ll^LlpR 



A3 -^"-^^ A3 



„HB _ HB 



^3 {WH)B^, X i^Ra^-'i^L + -j^{WH)B^, x i^L^^-'i^R 



„HW „HW 

^j:fT / ^'ttnt^.^" t ,,,, , 



^J^{H^a^H)W;, X i^Ra^^'^i^L + ^{H^a^H)Wl, x ^l<t^''^r , (11) 

where T) n ^ and Qf stand for each of the following gauge-invariant combinations of SM fields 

Of = V^H^V.H, {H^H)\ B^,B^\ W^^Wf , G;,G^^ , B^^B^^^ W^.W^. G^l^G^ ■ (12) 

But for the last three CP-odd combinations, the Qf monomials are precisely the dimension-four gauge and 
Higgs couplings of the SM Lagrangian. Operators with partial derivatives acting on -0 are systematically 




Figure 3: Scales A (in TeV) accessible using the rare K and B decays, for the vector current c]^^ and the scalar 
current cf ^. The plots on the right are obtained when MFV is imposed on the Wilson coefficients. 



discarded. 

The leading dimension-five operators involve either B^i^ or WH. The Bf^^ operators effectively assign a 
(mass-dependent) hypercharge to the invisible fermions, since under partial integration 

B^, X 4'B^<y''"i^L = 2iB^ X ^Pn^f^^L + 2m^B^ x (V^lT^V-l + i^Rl'^^n) ■ (13) 

With m^/A ^ 1 but otherwise arbitrary, the second term describes millicharged fermions, a scenario to be 
discussed later on. The H'^ H operators effectively correct the i}} mass after the electroweak symmetry breaking, 
with (setting to one) 

5m^ = wVA . (14) 

Taken at face value, naturality would thus prefer A > 24 (240) TeV if « 5ni^ <mB/'2. {jtikI'^-, as required 
to produce these states in rare B (K) decays. As seen in Tables [3] and SI this would push the NP effects from 
the FCNC operators ([5]) just beyond accessibility. Of course, these numbers are only indicative, a specific model 
need not produce these dimension-five operators, and even if generated, a sizeable deviation from naturality 
cannot be ruled out^ 

All the other Hf^ff. operators are, for our purpose, only marginally relevant. As explained in Sec. (TJ they 
reduce to the four-fermion operators of Eq. ([5]) once Higgs and gauge fields are coupled to quarks, with the 
four-fermion Wilson coefficients obeying the SM scalings ([7]). For example, the dimension-six operators with 
W V fj,H reduce, after the electroweak symmetry breaking, to an effective Z coupling to -0 which can be treated 
in perfect analogy to the SM Z coupling to (see Fig.[T]). The rare decays proceed through the flavor-changing 
hadronic Z penguin, and the four-fermion operators scale like c^'^^j^^^ ~ 9'^ll{rr)^^'^ with k^"^ given in Eq. ([7]). 

The final class of operators involves only one dark fermion field, and thus, from Lorentz invariance, violates 
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either lepton (C) or baryon (B) number: 

r^-C _ f>A£ _ A£ 

^AB,A£ = X ^rL + ^{H^H)H X V^flL + -j^B^^H x ^^a'^^L + x ^fla^'^L 

c^'C - - c^'-- - - C^^ - 

+ -^EL X ^PrL + -^DQ X {(jrL + -ttDD^ x ^jrU^ + h.c. . (15) 
A"^ A"^ 

For simplicity, the possible Dirac, SU{2)l, and flavor structures are understood, as well as the operators with 
V'-R ~^ V'l ■ Operators with three dark fermion fields are at least of dimension seven, hence not included. 

As explained in Sec. [TJ the phenomenology of the AB or AC operators is different from that induced by 
those of Hj^aj and 'H^'J . Since c"^® and c^^ cannot be simultaneously large as a tree-level ■0 exchange would 
induce proton decay, and since c^® cannot be probed with the rare FCNC-induced decays considered here, let 
us concentrate on the A£ operators. The renormalizable c^^ interaction, as well as c^^ after the electroweak 
symmetry breaking, mix ip with v^. As a result, these operators are bounded by neutrino masses, ipu behaves 
effectively as a right-handed neutrino, and the c^^ effective interactions include both FCNC and charged- 
current interactions. Since the renormalizable operator (the so-called neutrino portal) is not suppressed by the 
NP scale, its coefficient c^^ must be tiny. So, the effective interactions tuned by c^§' can be sizeable only if 
the NP dynamics responsible for such a suppression does not apply equally to all the A£ operators. 

In this respect, the semileptonic operator tuned by may be the most likely to escape such a suppression. 
It is the only one directly accessible with the quark FCNC transitions. However, the processes induced by the 
charged current component of the cf^ operator, J^(l— 75)u'^x-!/i(l— 75)^-'^, may give tighter constraints. Indeed, 
since there is no possible interference with the SM contributions, this operator enhances all the kinematically- 
allowed (semi-)leptonic decays of the K'^ , B+, and £)"•", as well as hadronic r decays. In particular, the most 
interesting channels are the K'^, £)+ — > transitions, whose rates are 

r«V^.>»."'-^'^--''-*''^'"-''-'*'f )7#(c?^H' . (16) 

647r \mai+m^jj \ ' ' J ' 

where P = B,K,D for I,J — (3, 1), (2, 1), (1, 2), = rrii/mp, fp is the decay constant (defined as 
in Ref. [H]), and the kinematical function A is defined in Eq. ([M]). Indeed, the SM charged-current Fermi 
operators are vectorial, hence these transitions are helicity-suppressed in the SM (proportional to m^isr). On 
the contrary, the relative strengths of the P~^ — ^ £~^tp transitions for £ = e, ^, t depend only on that between 
{cf'Y^^ , K —1,2, 3. For example, if (cf^Y'^^ is independent of K, the NP effects would be most easily seen 
for final states. 

Such an enhancement may be welcome in the B sector. The persistent tension between the Belle and BaBar 
measurements of B{B — )■ tv) = (1.68 ± 0.31) • 10~^ P^BHTO] and the predictions of the global CKM fits within 
the SM (see Ref. [T7] for a recent review) can be addressed provided < 3.5 GeV. Assuming (c^^)^^'^ ~ 1, 
the reconciliation of the discrepancy of the order AB{B — tX) ~ 10"'' would point towards A < 5 TeV, where 
the equality is reached for = 0. Note, though, that a non-universal lepton flavor structure in {c^'~'Y'^^ is 
necessary, in order to circumvent the current bounds from B ^ ev and B ^ fjLV aX the level of 10~^ [T^ . 

Actually, if the c^'' couplings are fully universal in their quark and lepton indices (i.e. {cf^Y'^^ ^0(1) 
for all I,J,K), and if ^ rriK, the best constraints currently come from the K£2 universality tesl[2l, where 
the NP effect gets enhanced by the small electron mass: 




rs*^(x^ez.,) + r^^(i^ 

From = 2.477(1) • 10"^ [H] and i?^^ = 2.487(13) • 10"^ we require - < 0.013 • 10"^ which 
translates as A > 82 TeV for ^ ttik- This is comparable to the scales probed with the FCNC modes. 

^The n(2 universality test is also constraining when < (we thank D. Bryman for pointing this out). Up to trivial 
substitutions, R,r is given by Eq. I I17I I even when tt — > ^tp is kinematically forbidden. A future measurement of R^^ at the 10~^ 
level [20] would give A > 67 TeV. Though lower than from Rk, the two are equally sensitive if (c^^)''"^/(c^^)'*"^ ~ Vus/Vud- 
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Some models with light weakly coupled fermions 

The operator basis ([5]) can be used to describe the SM transitions to neutrino pairs, 'ipL = Vf-^v^^Vr, 
t/jR = 0, by setting ah the c,; to zero but for cY,^. The SM rates for K decays are given in App. IA.3.1[ and those 
for B decays in App. IB. 3.11 With ipn — 0, aU but the vector FCNC operators c)^^ and c)^^ drop out. Since 
most NP models do not modify the particle spectrum below the electroweak scale, their effect simply enter as 
corrections to these two coefhcients pTll^ , and lead to the same vector-like spectrum as the SM contribution. 
Note that contrary to true dark fermions, the neutrinos are not neutral under SU{2)l ® U{1)y, and thus the 
basis ([H]) is not complete since it is lacking the charged current semileptonic operators. Said differently, if 
neutrinos had not yet been discovered, charge-current interactions would offer better windows. 

The presence of a light right-handed (sterile) neutrino corresponds to i/'i^ — vr, i/jl = in Eqs. © 
and (|15p . Of course, once is identified with i/^, it can be given a charge under U{l)c so that all these 
operators become A£ — 0. Some NP dynamics is nevertheless needed to couple vr, a gauge-singlet, to the 
quark currents. This can arise in the :/MSM model [53], based on the c^'-' coupling and with the quark flavor 
transition induced by the SM weak interaction. Alternatively, in some left-right symmetric models [24], the 
cf' operator can arise from combined Wl and Wr boxes. Note that these two examples are matched onto the 
operators of Eq. for which lepton universality tests are competitive with FCNC decays, as our analysis of 
the previous section shows. 

In the MSSM, the lightest neutralino xi is another electrically neutral weakly interacting particle which 
could be pair produced in rare decays. Though this particle is not neutral under the SM gauge group, it must 
be produced in pairs when R parity is conserved. Because the neutralino is a Majorana fermion obeying 

the corresponding operators disappear. In the general MSSM, there are then tree-level processes contributing 
to the combinations c)^^ -I- c^j^, c^i^ + Cjir, and to the scalar and pseudoscalar currents , X,Y = R. The 
flavor transitions are tuned by the off-diagonal entries of the down squark LL, and LR mass terms, and 
the overall scale A is set by the exchanged down squark mass. This (Class II) scenario was analyzed in detail 
in Rcf. 25], to which we refer for more information. 

Another possible new fermion is the axino a, the fcrmionic superpartner of the axion [26 . Depending on the 
model, it could be light enough to be produced in K and B decays. Note, though, that its Lagrangian couplings 
are flavor-blind, and further, involve the superpartners of the SM particles when R parity is conserved. So, the 
effective interactions are not only suppressed by the large scale A = fa, but also by the sparticle mass scale, 
loop factors, and the already tightly constrained flavor-violation occurring in the squark sector. We will not 
consider this scenario further because, with scales fa above 10^ TeV [52] (or even much higher in some models), 
the fermionic operators are far too suppressed, and signals should be more readily accessible using the single 
scalar axion production discussed in the next section. 

A flnal example is the dark sector millicharged fermion ip^ ^27]. Typically, these fermions arise when 
there is a new dark U{1) field coupled to the SM U{1)y through the kinetic mixing eB^'^V^^, as well as 
to some new fermion states initially neutral under the SM gauge group [55] ■ After diagonalizing the two U{1) 
fields, the dark sector fermions end up coupled to the SM photon, but with an arbitrary electric charge ee. 
Alternatively, this scenario could follow from the B^i, operator in Eq. p^ . 

When m^^ < me, various very tight cosmological and astrophysical bounds hold on e [57], ruling out any 
signal in meson decays. On the contrary, for m^^ ^ 0(100 — 1000 MeV), the bounds are much less tight, with 
e as large as 10~^ still possible. Since 4'e couples to quarks through the photon field, its coupling is fiavor-blind 
(Class III or IV of Table [5]) . The bounds in Table [3] are adapted to this scenario by setting A = M\y and 
rescaling the Wilson coefficients as c^"' — >■ ee^k^'' with k^"^ given in Eq. ([7]). Alternatively, the physics reach 
can be more accurately estimated by looking at the K and B radiative modes with a Dalitz pair, and rescaling 
their branching ratio by e^, 

B{P P'lpeipe) « X BiP P'l+r ) . (19) 

Up to simple phase-space corrections accounting for mi ^ m^^, this is valid as long as the Z boson does 
not dominate. From this, no B decay rate appears large enough to reach the interesting e < 10~^ region. 
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For example, with B{B {K,K*)e+£-) in the lO""^ range, a bound on B{B {K, K*)ij,ip,) at the 10~" 
level would be needed, far below the experimental prospects. Similarly, in the K sector, assuming about 10^"^ 
kaon decays, the only mode for which one could theoretically reach down to e w 10"'^ is — j-ipeipe, since 
B{Kl 7e+e-) = (9.4 ± 0.4) • 10"^^ pjj. 

So, the rare FCNC decays do not appear competitive when ip^ couples to ordinary matter exclusively 
through the photon. But, let us stress that if the couplings to matter of the SM and dark U{l)s are not 
perfectly aligned, the CKM suppression of Eq. ([7]) may be evaded, and rare decays would become prime sources 
of information. Also, as we will see in Sec. 14.31 in case the dark photon has a non-zero mass, it may be directly 
produced and competitive constraints can be derived. 



3 Invisible spin-0 boson 

If there is a light scalar particle neutral under the SM gauge group and under any dark gauge symmetry, it can 
be produced alone. The simplest effective operators are then of dimension five 

nLt = 4^Q7mQ X + ^D^,D X + ^hWQ x + ^HQD x . (20) 

No operator involving the tensor current or the SM field strengths arises, as these would require more derivatives 
and thus would be suppressed by additional factors of 0{mK,B /^)- Actually, the first two operators (as well 
as (Q pQ)(j), not explicitly included) collapse to the third and fourth upon using the tree-level quark EOM [I] 

i pQ = YlUH* + Y\DH , ipU^ YuQH*^ , ipD ^ Y^QH^ , (21) 
i.e., at the cost of the chirality flips 

^{4,Y'^ = (y,)"(4j" - {4^Y\Yay'\ ^i4^)" = (4^y^(Yt)" - i4,y'\Yiy'' . (22) 

Though in the rest of the paper, the EOM are always enforced, we prefer to keep all three operators here 
because they correspond to well-defined scenarios. On one hand, derivative couplings are characteristic of non- 
linearly realized symmetries, while on the other hand, the DQ x (f> operator is effectively a dimension-four 
Yukawa-like coupling after the electroweak symmetry breaking. 

If the scalar field is charged under some dark sector symmetry, or if the Z2 symmetry under — > —0 is 
imposed, it must be produced in pairs. The simplest effective operators are then of dimension six: 

'^tfat = ^<37^Q X iq^^^^^ + ^D-ff^D X z0t + ^hWQ X q^^q^ + ^HQD X . (23) 

^ y 

The minus sign of d in the first two operators is required as the plus combinations reduce to the third and 
fourth operators by partial integration and use of the quark EOM (PT|) . Tensor currents are not included since 
they are suppressed by a factor of 0{rnK,B / 

Bounds are derived on the operators involving quark currents of definite C and P, tuned by the couplings 

9s,p - ^ , gv.A ^ , gss.ps - ^ , 9vv,Av ^ , (24) 



where flavor indices are understoocQ. The corresponding differential rates are listed in Appendix IA.4I (jB.4p 
for K [B) decays. The physics reach is summarized in Tables [5] and |6] for one scalar in the final state, and 
in Tables |3] and S] for two scalars. As expected from Table [TJ the lower dimensionality of the operators (I^Ul) 
and translates as much higher accessible scales compared to spin 1/2 final states. 



*Note that the couphngs are hermitian matrices in flavor space (see Eq. lO), while (cg'^'*)''^ = (cg^^)^"' 
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Table 5: Production of a single invisible particle: Scales A (in TeV) accessible for the various operators, 
assuming bounds on the branching ratios of 10"^" for each mode. As for Table [31 the ranges of accessible 
invisible masses indicated in the first line are indicative, as the bounds are derived setting mx = (except for 
those channels which vanish, denoted by (*), for which mx = 100 MeV). These scales naively decrease when 
mx increases due to the phase-space suppressions, though the experimental acceptances need to be taken into 
account (see App. lA.ip . For the production of an invisible vector, see also Fig. |4l 
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Table 6: Production of a single invisible particle: Scales A (in TeV) accessible for the various operators with 
present (future) measurements (see App. IB.1|) . We assume mx = everywhere, except for those channels 
which vanish, denoted by (*), for which mx — 2 GeV. For the production of vector states, see also Fig. [5] 

The simplest effective operators involving Higgs and gauge fields are 

^' H^H X + X + j^H'^H X 03 ^ (25) 

00, H 00, G 40 

A' H^H X 0t0 + ^^H^ Vf^Hxcp^ d^'(j)+ ^^Q? x + ^H^H x (0+0)2 ^ (26) 

where Qf_^ are defined in Eq. (|T^ . Only operators of dimension less or equal to those in Eqs. ([20)1 and are 
considered. The possibility to write renormalizable couplings between and H embodies the so-called Higgs 
portal. One consequence is a mixing between and H, which has already been investigated in details, see e.g. 
Ref. [S], and will not be considered further here (but for a comment on A' in the next section). The other 
operators are a priori subleading compared to the Higgs portal, simply because of their higher dimensions. 
Note, though, that when is pseudoscalar, most of the Tif^'^ couplings drop out in the CP limit, leaving only 
those to Fp_i,F^^'^ with any one of the SM field strengths. 

Finally, the presence of a neutral scalar field does not open new possibilities compared to the SM to construct 
AB and AC couplings. The simplest operators are obtained by multiplying by (or by 0+0 if is not neutral) 
those of Ref. 0, and are either of dimension seven and AB = AC = 1, or dimension six and AB — 0, AC = 2: 

AC=2 _ AB=AC _ _ Ai3=A£ 

^Ae,A£ = -^HL^LH X + ^ ^3 Q^Q x L x + ' ^3 a,Q x Q^a'L x 

AB=A£ _ _ AB=A£ _ _ Ae=A£ 

+ ^-p^ — Q'^Q X U^E X -f — D^U X U^E x + — L x D'^U x + h.c. . (27) 



nj<P,3(t> 
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The A£ = 2 operator can produce an invisible ULVLcj) final state. However, not only is this operator of higher 
dimension compared to those of Eq. (PH)) . but its contribution to the rare decays proceeds through a neutral 
Higgs penguin (Class III in Tabled]), and is thus extremely suppressed. The AS — A£ = 1 operators have 
no impact on the AS = rare FCNC decays considered here. In addition, if < mp^n, those involving 
light flavors can induce AB = 1 proton or neutron decays. In that case, without highly non-generic flavor 
structures for the c^^^^^, the scale A must be close to the Planck scale. On the other hand, if > mp_„, 
these operators could still induce exotic B decays into an odd number of baryons plus missing energy. 

Some models with light weakly coupled (pseudo)scalars 

Typical non-linear symmetry realizations lead to derivative couplings of the type to the dark scalar 

field, with A given by the symmetry breaking scale F. Well-known examples of such light scalar states are the 
eLxion, resulting from the breaking of the PQ symmetry [F, or the familon, originating from the breaking 
of some family symmetry |29 . Only the latter naturally lead to FCNC couplings, since by design the family 
symmetry relates the three generations. For the axion, the dominant effect comes from its flavor-blind coupling 
to light quarks (Class IV in Table [2]). This is dominated by long-distance effects, speciflcally by the mixing of 
the axion with the light neutral mesons, as analyzed e.g. in Ref. [50] to which we refer for more details. Let 
us mention also that in some axion models, there is no direct coupling to light quarks, but rather couplings to 
the QED or QCD held strengths of the form (t>Ff,^F^"' or (pC^^G'^'' , as included in Eq. (^5)) . 

Similar derivative couplings arise from models of meta-stable supersymmetry breaking [3T], where 
a light pseudo Nambu-Goldstone boson (j) — P may be present. As discussed e.g. in Ref. [32], through the 
exchange of three gauge bosons, axion- like effective couplings to quarks are generated with Cy^ ~ ctfi^) and 
ai either the weak, strong, or hypercharge coupling. The scale A is the supersymmetry breaking scale, which can 
be around 10— 100 TcV in gauge-mediated supersymmetry breaking scenarios. These couplings are dominantly 
flavor blind, so according to our general rule (0, the flavor-changing vertices scal^ as Cyjfyj^ ^ a^{A)k^'^ . 
Thus, rare K and B decays seem unable to reach the interesting range A > 10 TeV. If supersymmetric particles 
are allowed to propagate in the loop(s), the operator may be enhanced, though the tight constraints on the 
flavor violation in the squark sector probably limit the accessible scales to A < 10 TeV range. 

Another scenario involving dark light scalar particles is the Next to Minimal Supersymmetric Stan- 
dard Model (NMSSM). But, in this case, the scalar mass is often larger than uib L33j, the scalar may decay 
too fast (e.g. to i^£~) to be considered as an asymptotic state in K and B decays, and the flavor transitions 
need to be induced by the (supersymmetrized) weak interaction. When MFV holds, the effective operators are 
then very suppressed (Class II scenario, see Eq. For all these reasons, the NMSSM does not appear as a 
likely scenario where our effective operators could play a role. Still, our formulas for the differential rates can 
be directly adapted to that case, and improve on the analysis of Ref. [M] by including more observables, and 
by a better treatment of the hadronic matrix elements. 

A more generic example is the singlet scalar model for dark matter, denoted S. In its simplest form, it 
enforces the Z2 symmetry and includes only the renormalizable coupling A' of Eq. (j26p , with = S* a real held. 
The capabilities of rare decays in probing this scenario were discussed in Ref. |35j . As described there, being 
flavor blind, the weak interactions are needed to induce the Cg^ coupling from the A' one (which couples 
SS to tt through a tree- level Higgs exchange). So, in this case, the key to interpret the numbers in Table [3] is 
the identification c^^^^/^/A^ — > X'k^-^ /m\, with k^'^ given in Eq. ([7]). See Ref. [3S] for more details. 

As a final example, the new invisible scalars could be the sgoldstinos S and P, the scalar superpartners 
of the goldstino [33 , for which the scale A corresponds to the fundamental scale of supersymmetry breaking. 
If light enough, these scalars could be produced alone or in pairs in K and/or B decays. In the former case, 
the coupling is of the form of c^^ g^, and is tuned by the LR squark mass insertions. In the latter case, the 
derivative interactions are of the form of Cy^ and Cy'^ with <p = S and ip^ — P, and are tuned by LL and/or 

^Note that this estimate differs from Ref. 1321 . where the flavor- violation is required to arise at the scale A from effective 
dimension seven operators. As explained in Sec. [T] dressing the effective dimension- five flavor blind q-f^qd/iP coupling with a W 
exchange allows for an electroweak-scale GIM breaking, so that the operator retains its 1/A suppression. 
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RR squark mass insertions. Note that since these squark mass insertions are rather tightly constrained by 
visible K and B observables, the accessible scales are limited by the MFV rescalings (O. 



4 Invisible spin-1 boson 

Compared to scalar and fermionic invisible particles, the presence of an invisible vector particle (denoted V) is 
significantly more difficult to parametrize. A consistent description of a neutral massive vector boson coupled 
to SM fermions is notoriously delicate. At the same time, adding a U{1) gauge group to the SM is one of the 
simplest possible extensions [28,37,38 . Furthermore, various theoretical models, whose motivations originate 
for example from strings |39) . extra dimensions |40j . or dark matter theories |41j . predict such new long or 
medium range forces, i.e. with masses in the MeV to a few GeV range. Such a dark vector boson would 
have many phenomenological implications, and has been intensely investigated recently |42] . It could also be 
produced in rare B and K decays, where it would show up as missing energy if sufficiently long-lived. In 
this respect, most models do also induce a coupling to leptons, but as long as the V — >■ vertex is not 

significantly larger than qq — >• V, our bounds should hold since producing a Dalitz pair through a virtual V 
exchange would push the rates well beyond the experimental reach. 

There are several ways to deal with light vector states, which we organize into three scenarios. For the 
first, the simplest FCNC operators involving the dark vector field are constructed. Those lead to decay rates 
diverging in the mv — > limit. This singularity is then treated assuming some kind of Higgs mechanism takes 
place in the dark sector (in a way similar to Ref. [38]). For the second scenario, the vector field is supposed to 
couple to SM matter fields only through its field strength, in a (dark) gauge-invariant way. This automatically 
ensures a safe mv — >■ limit, but significantly increases the dimensionality of the effective FCNC couplings. 
Finally, the third scenario considers only low-energy effective couplings of V to conserved flavor-blind quark 
currents [57] , as can arise for example from the couplings of V to SM gauge fields through the kinetic mixing [55] . 



4.1 Simplest FCNC operators 

For the first scenario, we consider the lowest-dimensional flavor-changing operators involving a vector fleld 
neutral under the SM gauge group, which are simply 

nlat^ = elQ^^Q xV^ + elDj^D x y^ (28) 

Thanks to the Lorentz condition dfj,V^ — 0, the leading correction starts at ©(l/A^). 

The flavor-changing quark currents Q^"i^iQ'^ and D^^^D'^ are not conserved when I ^ J (see Eq. (HI])). As 
a result, a naive rate computation with the polarization sum 

pol V 

diverges as my — ^ 0. This well known phenomenon is related to the impossibility of defining consistently 
the massless limit without an active gauge symmetry. It is interesting to compare this to the behavior of the 
dimension- four operators originating from the Z penguin in the SM (see Fig. [IJ. There, the Ward identity is 
violated only once the electroweak symmetry breaking takes place, and 



2 Q'l^Q' X H^V^H ^ " .^ 9v^Q'lt.Q' X . (30) 



{ATrAsM^^ '^^ (4^A 



SM) 



With Asm ~ v, this effective interaction is simply proportional to 5^, up to loop factors. However, the mass 
of the Z is also of 0(w), and thus it can never be on-shell once using the effective operator formalism. Instead, 
dimension-six four-fermion operators proportional to g^/M'^ ^ 9^ /v^ are relevant at low energy, see Eq. (jlj. 

In our case, V is light enough to be produced in meson decays. If it also gets its mass through some 
spontaneous symmetry breaking, theijfj my ^ s\ j^Vdark for some v^ark presumably similar or larger than v. 



® Actually, we should write my ~ ev^ark a-nd R ^ ^'^'^ some e and n > 0. The simplest situation corresponds to n = 1. 
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But as long as Vdark > 0, the limit my — > requires ^ — 0, which never diverges for physical processes. 
Actually, enforcing e)^^ ~ my/vdark in the my — )■ limit, any decay rate behaves as 

r(P ^ P'V) = [ d^p,vielj,MnV,Ak){elnMr / d<Pp,yM'^^M''* , (31) 

2mp J 2mp J vf^^^^ 

with Al^ the hadronic matrix element {P' \Q'' ^ ^iQ'' \P) or {P']^ ^^D'^P) and d<i>p'y the phase-space integra- 
tions. As expected from the equivalence theorem, this is precisely the rate one would derive from the axionic 
operators Q^jfj.Q^ x 9^<^ and D^j^D'^ x d^^cj) considered in the previous section, Eq. (pOj) with A = Vdark- 
For the pair-production of dark vectors, the simplest operators are of dimension six: 

r^V _ ^ VV _ VV _ VV 

+ -JT'Dj^ 'DuDx V^'V' + -f^Dj^D X V^V^" + -f^Dj^D x V^V'^" + -j^HQD x V^V^ , 

(32) 

with the field strength V^^ = - d^V^ and its dual V'"' = \e^'''P''Vpa. Note that those involving V^"" can be 
reduced to dv{Q"ffj.Q) x V^V and di,{Dj^D) x V^V, which are orthogonal to those tuned by c)^^ (which 
are missing in Ref. [6 ). These operators are given for completeness, but will not be considered further for two 
reasons. First, in the present minimal theoretical setting, there is no reason for the renormalizable operators of 
Eq. ((28)) to be absent, and those would clearly offer better windows. Second, the leading operators of "H^^JI] 
are those tuned by c^Y sb. since the others are comparatively suppressed by mK,B/v. But these operators also 
arise from the WDQ x V^^V^'^ and HQD x V^^^V^^^ operators considered in the next scenario (albeit in a 
rescaled form, see Eq. (|55| ). As explained there, these operators could become leading in the presence of a 
non-abelian dark gauge invariance, which would forbid single vector production. 

Phenomenology 

Let us consider separately the vector £v{d^ ^^j.d'') x and axial- vector EAid^Jfj.^dd"') x couplings. When 
my 7^ 0, both terms of the polarization sum (j29p contribute, so experimental bounds translate as constraints in 
the (my,£y) or (my, ea) planes. As explained above, we identify the regions where 'mv/£v,A > v as physical 
(for a similar reasoning, see e.g. Ref. [38]). 

As shown in Figs. S] and [SJ the pair (my,ey) is very constrained in both the K and B sectors, with typically 
mv/ev ^ V. Indeed, because the second term of Eq. ([^ is growing as l/rriy in the my — >■ limit, ey has 
to be correspondingly tiny to pass the bounds on the branching ratio. Said differently, the dark spontaneous 
symmetry breaking scale Vdark in Eq. pip has to be much larger than the electroweak scale. Interestingly, this 
remains true even when the flavor-violating part of £y satisfies MFV (green regions in Figs. |4]and[5]), i.e. is 
rescaled as ey — )■ jV^jVijIey, see Eq. ©. So, even if there is no theoretical requirement for the hidden symmetry 
breaking scale to be very different from that of the visible sector, current FCNC constraints nonetheless require 

Vdark > V. 

For comparison, the right panel in Fig. 0] shows the constraints one would get from a similar bound at the 
10~^° level on the K -fV channel. As this process is induced by the anomaly, the coupling is of the form 

eiV 

Cano^ = j^^{Re{ev)KL +ilm{ev)Ks)Ff,,V^''' , (33) 

with F^'^ the QED field strength. So, the current is effectively conserved, the l/rriy term in Eq. ([29]) drops 
out, and the my — > limit becomes smooth. 

4.2 Gauge invariant FCNC operators 

For the second scenario, we impose that only the dark field strength V^i, = 9^V^ — duVti and its dual V^^ = 
^e^'^'"^Vpcr occur in the flavor-changing couplings. This restores gauge invariance and ensures a smooth my — >■ 
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limit, but increases the dimension of the simplest operators. Depending on whether the vector is assumed 
abelian or non-abehan, the lowest-dimensional operators are either dimension six, 

fV _ J. JV _ - 

+ ^^7^ V.Dy^ + -^D^^^D x d^V^^^ + ^HQa^^D x V'" , (34) 

or dimension eight, 

vv _ vv _ vv _ vv 

^mlJlI] = -j^HWQ X y^,!/^"- + -j^HWQ X V^^V'^'' + -j^HQD x V^.V^-" + -j^HQD x V^.y'^'' 



VV .^v vv J. 

+ ^^^'^ ^-^^^ ^MP^"' + ^^7^ 25 X ^PP^"" + ^Djf^ ^ X ^PP^"" • (35) 

To reach this minimal basis, a number of identities and approximations were used. First, for the operators 

with a single vector field, the dual field strength V^'^ is absent from 'H^^^^[ll] since it can always be reduced to 

operators involving V'^'' using the quark EOM, integration by part, and the Chisholm identity |;1|. 

For the operators with two field strengths, we first note that summation over the generators of the adjoint 

representation of the dark gauge group may be needed to enforce gauge invariance. However, to keep things 

simple, we consider that after the dark sector symmetry breaking, only one (possibly complex) vector field 

is light enough to be produced at low energy, and simply discard the adjoint index. In addition, the non- 

abelian terms of the field strengths are systematically removed since the signatures we are after involve only 

two vectors. These restrictions, together with Vn.V'^'^ = —m^V^ and VnV'^'^ — 0, permit to reduce all the 

_ ^ — > 

operators with derivatives acting on the field strengths (as well as those involving Q^^{'D + V ^)Q). Finally, 
tensor operators involving two field strengths cannot be constructed since V^pVI is symmetric under fj, v, 
while H^D(jf"'Q X V^pVl reduces to H'^Daf^'^Q x V^pVP, upon using 2cr^'^75 = ie'^''"''^ a^fi. 

Reduction and phenomenology 

Not all the operators are relevant phenomenologically. First, we can discard all those involving a derivative 
acting on the quark fields, since they are comparatively suppressed by 0{mK,B /v)- Then, consider the operators 
involving d^V^^ in H^^jJII]. Upon enforcing the non-gauge invariant free EOM d^V^^^ = —myV^, they 
collapse to the dimension- four couplings of H^^JI], Eq. (1^5)) . with the identification 

^(Q7pQ, ^7p^) X a.y^'' ^ ^L^(Q^^Q, Dj^D) xV^ ^ -j^cZb^ ^ etj, . (36) 

This is thus an explicit realization of the H^^^p] operators, though with a major difference compared to the 
previous section. In the my limit, the coupling e)^^ here effectively scales at least like my instead of 
my, and the P — > P'V rates vanish when mv — > (compare with Eq. (PT|) ). The precise scaling is not fixed 
though, because that between c'Y ^ and mv is not known. If the vector gains its mass at the scale A from some 
spontaneous symmetry breaking, all we can say is that my ^ gA for some g, and c'^j^ ~ 5" for some n > 0, 
so that e)^fl ^ (my/A)^"*"". This scenario (with n — 1, i.e. g — ey — (my/A)'^) is shown by the light yellow 
regions in Figs. 2] and [5] Comparing these regions with those corresponding to the first scenario (in blue), the 
rescaling ([55]) is very expensive in terms of accessible scales. 

It is interesting to compare the reduction (j36p to that of the magnetic operator of HYnat [H] • After integrating 
by part and using the quark EOM (pij) . 

^HW'ap^Qj X V^-- = 2^H^H{{QYlY^pQ-' + D'^p{YlDy) x 

- 2i^hW^pQ'' X - 2^V''H^D'a^,Q' x . (37) 
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Figure 4: Above: Exclusion regions drawn from a bound on the n'^V (left) and Kl ^ jV (right) 

branching ratios at the 10""'^'^ level. The plots on the first line show the coupling 5 as a function of niv (in 
MeV) for the scenario I (blue, g = £y), scenario I with MFV (green, g = Vtdl), scenario II from the tensor 

operators sa^'^d x VJ^^ (yellow) and II from the vector operators 57^ d x d^.V^'^ (light yellow), and scenario III 
{g = ee). The grey area represents the region where my I g = Vdark < w ~ 246 GeV. The plots in the second 
line show the same, but replace g with A = mv/g (in TeV). 



The first two terms match those of H^j^jJI], Eq. (|28p after electroweak symmetry breaking. Since we started 
from a gauge-invariant operator, gauge invariance is now hidden in the quark-mass dependent relationships 
among the couplings e\ of 'H^^tP]; '^^H of those of the higher-dimensional operators in Eq. (j37p . It is 
only upon imposing these relationships that all the terms originating from the k^k^ /my piece of the polarization 
sum (j29p cancel out. Though this is another physically sound interpretation of the "H^^JI] operators, it is much 
easier phenomenologically to consider directly the operator H'^ a ^,^Q'^ x V^^'^ of "H^j^iJlI] from which they 
derive. 

Note, finally, that a similar reduction starting with the last two operators of 'H^^JII] can also be done. 
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Figure 5: Constraints on the coupling g against my from B — ^ KV (left) and B K*V (right) for the scenario 
I (blue, g = ev), scenario I with MFV (green, g = e|Vj^Vts|), scenario II from the tensor operators ba^^s x V^^ 
(yellow) and II from the vector operators 67^1 s x d^V^'^ (light yellow). The gray areas represent the regions 
where vny j g = Udarfe < f — 246 GeV. Lower plots: Same as above, but in terms of A = rny j g (in TeV). 



leading to the leading dimension-six non-gauge invariant operator of Eq. ((32)) 

yv _ 2 vv 

^H^DQ X l/^.y^"^ 3 -^Y^H^DQ x V^V" . (38) 

As explained before, in the present work, we consider this operator exclusively in the gauge invariant form, 
since the dimension-four couplings of Eq. (|28l) dominate if no dark gauge invariance is enforced. 

After the reductions described above, the only operators relevant for rare FCNC decays are those involving 
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the Higgs field, which we write after the electroweak symmetry breaking as (/ > J) 

+ '-j^d'd'' X y^.l/'^'' + l^d'-f^d-' X V^^V^'' + -j-d'a^^-f^d' X F^"^ + h.c. , (39) 

with, omitting flavor indices for simplicity: 

„ „v _|_ „v „, „vv I „vv „vv _|_ „vv 
«T.f = ^ 2 ' '^^^'P^ A 2 ' ^^P'PP A 2 ' ^ ^ 

The rates and differential rates for K (B) decays are in Appendix IA.5I (IB.Sp . while the physics reach are in 
Tables [3] to [HI Also, the constraints on the tensor operators are shown by the dark yellow regions in Figs. S] 
andini setting hrpf = {v/A)g and A = mv/g- 

4.3 Flavor-blind operators 

For the last scenario, the invisible vector boson is allowed to couple to conserved quark currents only, so that the 
singular term in the polarization sum (|29l) automatically cancels out in decay rate computations. To implement 
this, we have to drop the SU{2)l ® U{^)y gauge invariance requirement (see below), and couple V to quarks 
as 

= J^[c,] xy^ J^[c,]= J2 ^9^7m9- (41) 

q—u,d^s,c,b^t 

The Cq cannot be completely arbitrary but must reflect the flavor structures present in the SM. If that was 
not the case, then the FCNC couplings of the first scenario are in general present [33], and those would 
completely dominate in B and K decays. 

Using the MFV language, there are two possible fiavor structures. Either the Cq are universal or they are 
proportional to the Yukawa couplings. In the latter case, the top current would completely dominate, corre- 
sponding to Class III in the nomenclature of Table [51 In the former case, the remaining freedom corresponds 
to 

Cu = Cc ^ Ct = Cu , Cd ^ Cs ^ Cb = CD , (42) 

with cjj and a priori arbitrary. Note that this two-parameter freedom means that it is always possible 
to write J/i[cq] in terms of the electromagnetic and the baryon number currents only, as was pointed out in 
Ref. [33] • When the invisible vector is aligned with the photon, rather tight constraints are set by flavor-blind 
observables, e.g. the muon g — 2, quarkonium decays, beam dump experiments (see e.g. Ref. [3] for a recent 
analysis). While we will compare our constraints with those, let us stress that consistency does not require VJ^ 
to couple to lepton^. If it is leptophobic, many of these limits can be evaded. Hence, only the purely hadronic 
production oiV in K or B decays will be considered here, and should be compared e.g. with those from tt*^ or 
quarkonium decays with missing energy. 

As explained in Sec. [2 it is very different to probe the couplings to heavy (Class III) or to light quarks 
(Class IV), so we will analyze each situation in turn in the next two sections. But before that, let us return to 
the issue of an SU{2)l U{1)y gauge invariant origin for the conserved currents in Eq. (liTj) . 



Flavor-blind SU (2) l (E) U {1)y invariant couplings 

In Sec. 14. 11 the 1/my singularity of the decay rates arising from the FCNC couplings of Eq. ([55|l was interpreted 
in terms of a dark symmetry breaking scale Vdark ^ v. Specifically, we enforced that the coupling constant e of 
to quarks and its mass satisfy my ~ VdarkS, so that it is always the finite combination /my — >■ I/w^q^^, which 
occurs in decay rates. But, generic vector or axial-vector quark currents are conserved for massless quarks, 

''Anomaly cancellation may require the leptons to be charged under the dark gauge symmetry, depending on the detailed 
particle content and dynamics of the dark sector. We do not consider such constraints. 
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which is the case when SU{2) l'SiU{1)y is exact (for simphcity, we disregard the chiral current anomahes, which 



do not concern us here). So, the presence of the apparent l/rriy singularity could be due to the electroweak 



symmetry breaking instead of to the dark symmetry breaking. 

To make this statement more concrete (see also the discussion in Ref. [46]), let us write down the flavor-blind 
renormalizable couplings between V and SM gauge, Higgs, and quark fields (the flavor-blind summation over 
/ = f , 2, 3 is understood) 

nlJlll] = £iB^. X V^"" + e2iH^'^^H •KV^'+ sgB,,, x W'' + eniH^H) x V^V , (43a) 
^^^,,,[111] = eB{Q'l^.Q' + D'^^D' + U'-,^,U') x V + eoiD'^^D') xV^+ eu{U'j^,U') x V . (43b) 

These dimension- four couplings presumably dominate over the higher dimensional flavor-blind operators, which 
we do not list here (See Ref. [H]). 

Only the £i, £2, £e, and Eb couplings are compatible with a dark sector gauge invariance associated to V, 
since they vanish under — > (see Eq. (PT|) ) upon partial integration, and using the free Higgs boson and 
quark EOM. The eb coupling involves the (conserved) baryon number current, and is thus directly matched 
onto Eq. (|4fT) . The other couplings, Eh, £_d, and Eu, would break the dark gauge invariance and are thus 
discarded. Note that eh must in any case be tiny if V is to be light enough to be produced in rare decays. 
Indeed, the {WH) x V^V^ coupling gives a mass to V after the electroweak symmetry breaking, so barring a 
large cancellation between the dark and visible Higgs sectors, Eh S ™k(_b)/''^^ ~ 10^^(10"'*). We can further 
discard the Eg term, which is a total derivative and is relevant only for magnetic monopoles |47) . 

The two remaining couplings are £1 and £2. The first one is the celebrated kinetic mixing |28] . Since 
Bfj^ = cos Ow^ti — sin^vF-^/i, we can rewrite it as 

B^,^ xV"" ^2 cos 9wr^"' xV + 2 sin OwZ^ x df.V'"' . (44) 

The piece proportional to the (conserved) electromagnetic current d'^F^i, — — J^™ can be matched onto Eq. (|H|), 
while the other one mixes Z and V. The £2 coupling also generates a direct mixing x F'' after the electroweak 
symmetry breaking, since WV^H gv^Z^. 

The Z — V mixings induced by £1 and £2 are very different. The kinetic mixing Z^ x d^V^^^ is safe in the 
mv — ^ limit, because it is insensitive to the electroweak symmetry breaking. Actually, once the Z is integrated 
out, this Z^ X d/^iV^" vertex together with the SM flavor-changing hadronic vertex of Eq. ([50)1 generate the 
gauge-invariant operators c^^ of Eq. (IMl) . with c'^ ^ Eig^V^jVij and A ^ Mz- On the other hand, doing the 
same with the £2 coupling generates the dangerous operators of the first scenario, Eq. (j28p . with 

isD" ^ £2..^ X ^r^.VuVu - X -£,V:,Vu . (45) 

However, in parallel to the above FCNC operator, £2 also corrects the V mass as Sm^ ~ s^V^ ■ As a result, 
TTiy'^* _j. Q j-gquij-gg £2 — > 0, ensuring again the safety of all decay rates in the massless limit. This shows that 
indeed, the visible symmetry breaking scale can play the same role as a dark symmetry breaking scale. The 
only difference, besides v 7^ Vdark in general, is that the former relies on the SM dynamics to drive the FCNC, 
and thus brings in the loop and CKM suppression factors, see Eq. (|45)) . The allowed range of £2 values can 
be derived from the green areas in Figs. 2] and [SI up to the rescaling by g/lQn^ ^ 10"^. Note that £2 values 
acceptable for rare decays ensure that Sruy ~ e^v^ can be safely neglected. 

The above electroweak mechanism ensuring a safe massless limit may render the extension to a two Higgs 
doublet model desirable. Indeed, there would then be two different Eh couplings, and an additional conserved 
current can be constructed [Sj, whose charges are aligned with those of the Peccei-Quinn (PQ) symmetry [48] . 
Combined with the conserved lepton number current, this allows in principle to make V completely leptopho- 
bic |49| . The cost being the presence of a flavor-blind axial- vector quark current, not conserved at low energy. 
The corresponding 1 / my singularity is nevertheless under control thanks to the additional sources oi Z — V 
mixing present in the PQ current, i.e. tuned by the same coupling constant. We will not further elaborate on 
this construction, but just retain that a leptophobic setting is in principle possible, allowing to evade the many 
low-energy constraints based on the 67^6 x and /i7^A* ^ couplings ^4u42j . 
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Apart from the Z—V mixing effects, matched onto Eq. (1^51) or Eq. ([M)) . the Eb and e\ coupUngs are genuine 
SU{2)l (8) U{1)y invariant reahzation of the two conserved quark currents of Eq. (|1T|) . Let us now see how to 
constrain them from rare decays. 

Phenomenological constraints on the couplings to heavy quarks 

Once the heavy quarks are integrated out along with the weak bosons, the presence of V m K and B physics 
is felt through the operators of the second scenario, in particular '^^^^([11]. For example, the last operator in 
Eq. (|34p is induced in complete analogy to the electromagnetic operators describing 6 — > S7 and s — > c?7 in the 
SM. This situation is thus simple to account by adapting the coupling of 'H^^jJII] according to Eq. ([7]), and 
setting the scale A at Myy- Alternatively, a more precise estimate can be obtained when the new invisible vector 
boson is very light and aligned with the photon (in the quark sector). If we set cjj = 2/3ee and cd = — ee/3, 
the branching ratios for b sV and s — > dV are obtained by rescaling by £^ the SM predictions for the b ^ 
and s ^ d'j processes, up to simple phase-space corrections. 

Specifically, in the B sector, the branching ratio for b — > sV is 

B{b ^ sV) ^ \e\^ Bib ^ s-ff^ , B{b->sjf^ = {3.15± 0.23)- IQ-'^ [5D] , (46) 

when my <C mB and E^^y > 1.6 GeV. This cut on the photon energy is actually at the opposite end of 
phase-space compared to those set for b — ^ sX. But even without a definite prediction, it is clear that the 
expected sensitivity of about 10~^ in the B — ^ [K, K*)X channels would at best probe e down to a few percent. 
For comparison, typical bounds on e derived from flavor-blind hadronic observables are currently down to the 
10-3 range 

The situation is worse in the K sector, where only CP-violating observables are sensitive to the short- 
distance (c and t) magnetic operator. As analyzed in Ref. |51j . those are beyond experimental reach even in 
the SM case, and thus cannot be used to set constraints on e. This is actually clear from Table [SJ rescaling by 

~ IQ-^, the scale A ends up well below the electroweak scale. 

So, rare K and B decays are rather ineffective at constraining the presence of a new flavor-blind vector 
coupled exclusively to heavy quarks. Fortunately, in many cases, as e.g. from Eq. (j44p . universality holds and 
this vector must also couple to light quarks, where the situation is much better. 

Phenomenological constraints on the couplings to light quarks 

In this case, the CKM factors strongly favor the K sector to derive competitive bounds. At the K mass scale, 
only the u, d, and s quarks are active quark degrees of freedom. Adopting a matrix notation in the q = (?i, d, s) 
flavor space, "^^^^[111] takes the form 

nlfilll] = e g-7,Q'g x y^ Q' ^ eQ + e'l , e ^ e' ^ 2E±^ , (47) 

with 1 — diag(l,l,l), Q = diag(2/3, — 1/3, — 1/3), and e the QED coupling constant. So, from the point of 
view of low energy physics, there are only two possibilities: either is effectively aligned with the photon (e 
term) or its charges are proportional to baryon number (e' term) |44) . This Ti^jf [HI] coupling must be directly 
embedded within Chiral Perturbation Theory (ChPT) [S^- At the leading order, the field enters only 
through the covariant derivative acting on the meson fields 

D^,U = - leA^ [Q, U] ~ iV^ [Q', U] = d^,U - ie{A^ + eV^,) [Q, U] . (48) 

The e' term cancels out in the commutator, leaving coupled exactly like the photon A^. This ensures the 
absence of a direct K nV coupling at leading order, relegating them to 0{p^). Such a direct leading order 
coupling only exists when the d and s charges are different. Indeed, in that case, the generator Q' would no 
longer commute with that of the weak interaction. This is another way to see that when the universality (1421) 
fails, the dimension-four FCNC couplings of T-L^^f. [I] should be allowed. 
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Experiment (for = 7) [12] Indicative reach 
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Table 7: Indicative experimental reach of the radiative decays for a new light vector boson aligned with the 
photon (£) or with baryon number (£'), assuming about 10^^ kaon decays are observed and bounds in the 10"^^ 
range are set. For the K'^ — ^ tt^V mode, the reach in |£| drops as my decreases (with limits at the 10"^ level 
for mv « 50 MeV, 10"^ for my w 5 MeV, and only lO'^ for my « 0.5 MeV). 



To get bounds on e is rather immediate since the phenomenology is completely analogous to that of the 
radiative K decays (see Ref. [51] for a recent review) . It sufHces to consider the dominant radiative modes and 
replace one photon by V. When it is massless, the rates are obtained from those in QED simply by rescaling by 
(or by a' /a, if one defines a' = e'^a). When massive, the amplitudes are essentially the same, the polarization 
sum is identical (since the QED Ward identity holds), so the main impact is a reduced sensitivity due to the 
truncated phase-space. Assuming that about 10^^ kaon decays will be analyzed in the next generation of 
experiments, and bounds in the lO'^^ range are set, the reach in e is thus naively 



_2 ^ B{K -^mr + mj + V) 10" 



~ B(/s:^n7r + (m+l)7) B{K ^ mr + {m + ' ^ ' 

for massless (or very light) vector boson V ^ and n + m > 0, n < 4. Since the current bounds on e are down to 
the lO""^ range, competitive bounds could be obtained from all the modes with B{K n7r+ (m+ 1)7)) > 10^^. 
Those are listed in Table [7| 

The K — > ttV channel is special because K — ^ 777 is forbidden due to gauge invariance. Further, even when 
off-shell, K -> 7r7* vanishes at leading order in ChPT, and so does K — > irV . The leading contribution thus 
starts at 0{p^), from loops and local counterterms, and is approximately given by [53| 

A{K+{P) ^ 7r+y(g)) = £^a+ {q^P^^ - q^P ■ q) e;{q) , (50) 

with a+ an 0(1) constant. The Ks rate is expressed similarly in terms of the 0(1) constant as, while that for 
K2 ~ Kl is CP-violating and thus driven by heavy quarks (since lm{V*i.Vud) = 0). From this amplitude, we 
get the rate 



r(if+ ^ .^V) = ^^«lg^l«+ r ^ :j"'^ r'v , (51) 



with A(l, rj,ry) the kinematical function defined in Eq. (j95p . and = mi/mji- Normalizing with the K^ — > 
TT+e^e- process to get rid of a+, 

p2 Q\3/2(i 2 2 \ 

BiK+ ^ 7:+V) = - 4 xB(K+^ n+e+e-) , (52) 

with $e ~ 0.145 the K+ Tr+e+e" phase-space factor, and B {K+ Tr+e+e") = (3.00 ± 0.09) • lO"'^ [H]. 
Reminiscent oi K ^ nj, the rate vanishes in the my — limit. This seriously hampers the reach in |£|, as 
indicated in Table [7| Even in the most favorable window m^r < my < 2m7r, K^ ^ tt+V is less sensitive to |£| 
than Kl — > by about an order of magnitude, see the red regions in Fig. |31 

To get bounds on e' is more difficult because it cancels out from the 0{p^) Lagrangian, and thus also from 
the 0{p'^) meson loops (see Eq. (HHj)). It can thus occur only in some local counterterms involving the V^^ 
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field strength (of vanishing anomalous dimension since there are no divergent loop contributions), and in the 
odd-parity anomalous 0(j>'^) Lagrangian. The former are very suppressed compared to the loop contributions 
induced by QED and by the e piece, and will be neglected [SI] ■ 

Concentrating on the odd-parity sector, only the Kl — >■ jV mode appears useful to constrain e' (its 
anomalous amplitude, driven entirely by the up quark [54], is sensitive to both e and e'). The magnetic direct 
emission amplitudes in K ^ ttttj are significantly more suppressed and difficult to access experimentally. 
Rescaling Kl — J> 77 according to Eq. (HHj) shows that a bound on Kl — >■ at the 10~^^ level would probe 
couplings down to at least \e,e'\ < lO"**. Interestingly, this is more than an order of magnitude better than 
using the flavor-blind transition tt'^ — > •^V, for which the best limit is 3.3 • 10^^, i.e. \e,e'\ < 6 • 10^'^ [T^ . 
Further, the range of accessible V masses is evidently larger in K decays. 



4.4 Baryon and lepton number violating operators 

As for the dark scalar, Lorentz invariance requires an even number of SM fermion fields. However, the vector 
field index allows for alternative chiral structures compared to the Weinberg operators [2] . Specifically, keeping 
only operators of leading dimensions, 

•hXS = ^HL^V^LH X + ^HL^LV^H x ^HL^ LH x y^F" 

-f -^HL^a^,^LH X V^"' + h.c. , (53a) 



^AS=^A£ = -jfD b^-^D X ^7"^ X + -f^D'-b^^D X Lj-'Q x + h.c. , (53b) 
,v,vv 



A3 ' ■ ' A3 

„I „l J I 



^ab:ac = -fiJl^ X V, + -j^Ja X V,V^ + ^ Jr X V + h.c. , (53c) 

where the tensor stands for g^^ or ct^^, the set a = QQQL,QQUE, DUUE, DUQL denotes the gauge- 
singlet combinations of fields of the Weinberg operators , the corresponding currents are defined as 

Jabcd = A^B X C^'D , (54a) 
J^^P^ = A^VB X C^D, A^'B X C^VD, A^V^'B x C^D , (54b) 
Jabcd ^ f^^'' B x C^D, A^B x C^a'^'D, A^-i^'B x C^^D , (54c) 

and where the SU{2)l triplet contraction for the QQQL current is understood. Note also that c^ is anti- 
symmetric in flavor space, c\ is symmetric, while c\ 2 have no particular symmetry, though one of them is 
redundant when flavor-diagonal. For the tensor current, which of the three Dirac structures does exist depends 
on the chiralities of the fermions involved (and may require reordering the fields using Fierz identities). The 
only other possible vector current is J'^cjjll^ vanishes upon Fierzing due to the antisymmetric SU{2)l 

contraction of the two lepton doublets and thus requires two more Higgs fields. The superscripts I and II 
refer to the scenarios discussed previously, i.e. separates those operators which explicitly break a dark gauge 
invariance associated with V from those which do not. 

All these interactions have high dimensions, especially compared to the renormalizable couplings to SM 
particles of Eqs. (|28p and Phenomenologically, they may be relevant only when A is not too large, which 
requires the absence of direct FCNC couplings with V, see Figs. S] and [SJ At the same time, the AS = ± A£ 
interactions can induce nucleon decay when my < mrin (note that AS = — A>C interactions contain at least two 
d quarks and do not contribute to p'^ decay at the leading electroweak order), and thus require either the scale 
A to be extremely high, or the Wilson coefficients to have highly non-generic fiavor structures [TT| . 

A scenario with ms > ray > rUn is interesting since astrophysical, leptonic, and nucleon decay bounds are 
essentially circumvented. In that case, rare B decays into an odd number of baryons (together with an invisible 
V) may offer the best windows for the AB — ±A£ interactions. Note, though, that if these interactions occur 
concurrently to those of Eq. (pS]). (p4)) . or (|43|. the V may occur as an intermediate state, bringing back the 
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tight proton decay constraints. It remains to be seen whether in that case, signals in B decays are nevertheless 
possible. Such virtual exchanges are beyond our scope, since in the present work, we require the dark particle 
to be sufficiently long-lived to escape as missing energy in rare decays. 

Finally, as for the dark scalar scenario, the AC = 2 interactions can produce three-body invisible i/LiyLV 
final states, but are not particularly interesting for FCNC decays. Indeed, they are unable to induce quark 
flavor transitions, and the FCNC decays would thus proceed through an extremely suppressed hadronic Higgs 
penguin. 



5 Invisible spin-3/2 fermion 

Spin-3/2 particles are described by Rarita-Schwinger fields, denoted ^'p, which transform as spinors with a 
vector index under the Lorentz group. The corresponding Lagrangian kinetic term can be written as |55j 



In addition, these fields are also subject to the conditions — (spin-3/2 projection), [il/) — m<^)'^^ = (Dirac 
equation), and d^j,'^^ = (Lorenz condition). For external states, the spin summation is performed as |55j 

n(p)M. = E ^(p)^^(p)^ = + - \p^.pP..Yln , Pep = gap - (56) 

spin. ^ 

The sum over spin for v{pY^ spinors is given by — n(— p)^i,. 

We distinguish the possible operators for the pair-production of these fields by their Lorenz structures, 
which can be scalar, vector, or tensor-like (even though Eq. ((55)) is written down for a Majorana field, ^E* will be 
taken as complex from now on). Taking into account the above-stated conditions reduces the possible leading 
operators to 

= +^Q^pQ X ^'r^p + ^Dj^^D X W'r^p + ^Qj^Q X ^''7^75*p + ^Dj,D x ^''7^75*p 

S P S P 

+ ^HWQ X W^*^ + ^HWQ X ^'^75*^ + ^HQD x W^^-^ + ^HQD x ¥^75*^ 

T T T 

+ ^H^Da^.Q X ^I^*"-! -I- ^H^a^^Q x ^[^75*^! -|- ^nWa^^Q x ^.a^"-*" 

+ ^HQa^^D X + ^HQa^^D x *[^75*-l + ^HQa^.D x ^^a'^''*'' , (57) 

where I'ITvE''"! = i(*''r^''' - l'TvI''^)/2. Only the leading operators of each kind are kept; operators 
with additional derivatives are systematically discarded. For the vectorial couplings, the operators involv- 
ing e^'^'^'^'i'jyjp'^cr and e^'^'^'^'^i^-f^-fp'^^a have been reduced to the others using the Chisholm identity. Finally, 
tensor structures are similarly reduced using the Chisholm identity together with a^^e^^P"^ = —2ia'^'^^5, which 
in particular permits to get rid of the 6^^°''''^ a-Js'^ p and ^'''cr'^'^75^p structures. 

The effective couplings of dimensions up to six involving gauge or Higgs fields are easy to construct from 
those in Eq. (ITT|) . 

- +'fH^H X + ^-fH^H X ^75*, 

S P T 

+ ^Bp, X ^[^^^^l + ^B^, X W[^75*''l + ^B^, X ^pa^'^*" 

C"^ 

+ -^Bp, X ^[^^-''l + -^B^, X W[^75*''l + -^S^, X '^pa^"'^P 

Y A 

+ 2^^H^^^H X ^"7^*^ + gzijt Vpi? X ^''7''75*p . (58) 
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After the electroweak symmetry breaking, the situation is similar as for spin 1/2 fields. The dimension-five 
operators in the first line generate a correction to the 5* mass (upon enforcing the ^ = constraint), those 
in the second, third, and fourth line couple to the photon and to the Z boson. Note, though, that \1/ does 
not become millicharged in the usual sense, as these effective operators do not match those derived from the 
minimal substitution principle (which, in any case, is not consistent for spin 3/2 particles [56]). 

The last class is made of operators involving a single '3/ field. As for the spin-1/2 case, Lorentz invariance 
requires an odd number of SM fermion fields, so these operators break either baryon or lepton number: 

A£ _ A£ _ A£ _ 

^AS,A£ = ^^M^ X + X *'''7''li + -Jt'^U^^' X ^''^^''li 

„A£ _ _ A£ _ _ AB _ _ 

+ ^EcT^,,L X *['^7'^1l + ^Do^,,Q x ^^i^^^^L + —Do^.D'^ x ^^i^^^V^ + h.c. . (59) 

Notably, compared to the spin-1/2 case (ITSt . no renormalizable coupling can be constructed, and thanks to 
the extra derivative in the dimension-five operator, there is (of course) no direct mixing between ^ and vj^ 
after the electroweak symmetry breaking. Phenomenologically, the signatures for the AS operator would again 
require specific searches in B decays, while those for the A£ operators are to be found in semileptonic decays. 
Note, however, that the interactions in Eq. ((5^)) are more difficult to access than those for spin 1/2 invisible 
particles, Eq. because the tensor matrix elements vanish, (^\d}- a^'^ \P^) = 0. This means that c^'~ 

does not contribute to the P+ ^+5' decays, but only enters in the P — > P'l'^ decays for which the helicity- 
allowed SM contribution P — >■ P' (.v is large. Purely leptonic processes are only sensitive to higher-dimensional 
operators involving covariant derivatives acting on the quark or lepton fields, for example D'Dfj^Q x "^^L. So, 
the AiZ operators do not appear promising in searching for dark spin 3/2 particles, and will not be further 
considered here. 



Reduction and phenomenology 

The non-conserved quark flavor- changing neutral currents break the gauge symmetry appearing in the La- 
grangian (|55p when mif — )■ 0. As a result, the 1/mip terms in the spin sum (I56p are not projected out, the 
massless limit is singular, and we can force A up to arbitrarily high values simply by decreasing m^. The 
situation is analogous to that encountered for the massive vector case in Sec. 14. 1[ and may resolve itself in a 
fully dynamical theory in a similar way. To get physically meaningful bounds on the scale A, there are two 
possible routes. 

The first procedure is inspired from the supergravity setting |57j . where the spin 3/2 gravitino mass is 
related to the supersymmetry breaking scale as h^susY — iV3m<s/Mpiancky^'^ with Mpianck — (STrG'Ar)"^/^. In 
some sense, this can be understood as the fermionic equivalent of the constraint my ~ gvdark enforced for the 
vector bosons (which would here be insufficient given the harder (m^^)^ singularities occurring when there are 
two spin-3/2 particles in the final state). Indeed, when Asusy ^ Mpianck, "t-* is very small and only those 
terms originating from the singularity of the spin sum (1561) are relevant |55] , 

UipU "="^° ■ (60) 

This projects out the ±3/2 helicity states, leaving the ±1/2 goldstino helicity states in a way similar to Eq. (1511) 
for the massive vector boson. Specifically, the spin-3/2 operators become equivalent, in the — limit, to 
the spin-1/2 derivative operators obtained by replacing 

f2d^^ 



= V^— • (61) 

V 3 

Given that the effective operators are at least of dimension six, there are at least two powers of A = Mpianck 
which can be eaten away by enforcing Amip — >■ A^gugy- 
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The supergravity scenario is thus characterized by the rescahng Asusy = ('\/3?7i*Afpiancfc)^^^- So, even if 
4* is here not necessarily identified with a light gravitino [S^, let us assume that 

A ^ A^/mq, , (62) 

where A may not be related to Asusy in any way but could denote some dark sector symmetry-breaking scale. 
Phenomenologically, this rescaling permits to derive sensible bounds on A from the rare P — > decays, 
even when mif <C mp. 

A second route would start from a basis made entirely of gauge-invariant operators, i.e. involving the 
field-strength ^E*^^ = d^'ii, — and its dual XP^jy = e^^pad'''^'' . Though the rriif — > limit would always be 

smooth, we do not perform this construction explicitly because with two field strengths, there are too many 
operators for the basis to be useful phenomenologically. Instead, we simply remark that starting from such a 
basis of gauge-invariant operators, it must be possible to generate the 'H**^ operators by partial integration 
and use of the EOM, exactly as for the massive vector boson. 

Specifically, when only the spin-3/2 EOM is used, i'^f^'i?'^'^ ~ m^s/'i?'^ and 7p\t'^'' ~ —777,^75^'^, an extra 
factor m^/A^ is generated. For example, 

l^H^DYYQ X ^apl^r^^.u - "^H^DYl^Q x ^p*. , (63a) 

YiQlpQ X ^''"^''^7A.7^7"7^*/^. - ^Qlo^Q X *p7"*^ • (63b) 

This is analogous to the reduction p6p for the vector boson. By contrast, whenever the reduction involves the 
quark field EOM (PTI) . the gauge invariance ends up hidden in relationships among the Ci of exactly like 

in Eq. ([37]). The resulting operators are then suppressed either by the light quark masses, or by derivatives 
acting on the quark fields. 

Enforcing some cancellations among the operators is incompatible with our procedure of turning on one 
operator at a time. So, we consider only the situation where the P — >■ P'^f^* decay rates are finite in the 
— !■ limit thanks to the rescaling 

c^ ^ , (64) 

with Ci ^ 0{1). Note that this cures the singularity (pOj) . Away from the strict = limit, the other terms 
of the spin sum (|56p also contribute and tend to suppress the rates. 

To derive the bounds on the scale A from the rare decays, we must impose one of the above two prescrip- 
tions (I62[) or (j64p to make sense of the — > singularities. Comparing them, these rescalings appear precisely 
equivalent for the dimension-six operators of H**^. For the dimension-seven operators, the gravitino- like rescal- 
ing (|62p leads to an additional suppression by m^/A compared to (|64p. making them completely irrelevant 
(remember that is assumed smaller than tok^b)- Thus, the bounds we quote for these dimension-seven 
operators are understood to hold only for the second scenario. 

As for the other types of invisible particles, we rewrite the various operators in terms of currents of definite 
C and P by introducing the fourteen complex couplings (for each s — d, 6 — > s, and b ^ d operators) 

fxv,xA = (X = V,A), hs.xP = J-^^ (X = 5,P), = - £k±ck ^ s,P,T) . 

(65) 

The rates and differential rates are in App. IA.6I (IB.6P for K (B) decays, and the corresponding bounds are 
shown in Tables |3] and SI As these numbers show, the rescaling prescription pushes the dimensionality of the 
operators to eight or nine. For such dimensions, the accessible scales A are at or even below the electroweak 
scale, as expected from Eq. ([T]) and Table[Tl and the viability of the SU {2)l(^U{1)y effective operator formalism 
becomes questionable. Said differently, if the scale A is above the electroweak scale, the presence of a dark spin 
3/2 fermion should have no impact on rare FCNC decays. 
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Hmat Dim 


Hint Dim 


I-Lab.ac Dim 


"'"/^ tptp 




il) 6 


(HH) 5 


HSI 6 (4) 

(>8) 


^ 




6 


(123 3 
(PS)1 4 


^ — 8 (7) 


Direct, V 

1 

Gauge, V 
VV 




HMD 4 
(EH) 6 

(El 6 

(ESD 8 


4 
6 


7(7) 

m 

10 (9) 






dSTD 6 


dsg) 5 


dSSl) 6 (5) 

(> 8) 



Table 8: References in the text for the various pieces of the basis of effective operators. Dimensions are denoted 
with a bar when the leading operator involves a Higgs field reducible to its vacuum expectation value after the 
electroweak symmetry breaking. In the last column, the dimensions are indicated for AS operators, irrespective 
of their A>C components, while the dimensions in parenthesis are those of the purely A£ operators. For dark 
vector fields, we distinguish between direct couplings, for which the my — > limit is formally divergent, from 
those where a dark gauge invariance effectively survives (even though full invariance is not imposed, since we 
allow for mv > 0). Finally, the dimensions of the operators involving spin 3/2 fields increase when ensuring a 
sensible — > hmit, see Eqs. (|52|) and (|M)) . 

6 Conclusions 

In this paper, we presented a complete basis of SU{i)c ® SU{2)l ^ U{1)y invariant operators involving SM 
fields together with a yet undiscovered light invisible spin 0, 1/2, 1, or 3/2 state, neutral under the SM gauge 
group. As summarized in Table |H1 the operators are organized into three classes: couplings to SM fermions, 
couplings to SM gauge and/or Higgs fields, and baryon/lepton number violating couplings. We retained the 
operators of lowest dimensions separately for each class. As a result, most of them do not strictly qualify as 
portals since they are suppressed by the NP scale A. However, it makes sense to extend this denomination to 
those operators for which the experimental constraints push A far above the electroweak scale. For example, 
the typical scale for a dimension- five FCNC operator is greater than 10000 TeV (see Table [1]), while it can 
even be close to the Planck scale for those inducing proton decay. For this reason, in the present paper, we 
systematically investigated the FCNC operators, and derived bounds from the rare FCNC transitions. 

Our results can be split into two parts: those concerning the basis of operators itself, and those related to 
its phenomenological impact on the rare decays. Starting with the basis, some of the main features are: 

1. First and foremost, it must be stressed that even though we concentrated on the rare K and B decays 
involving missing energy to derive bounds on the operators, our basis is completely general and could be 
used equally well to investigate signals e.g. in lepton flavor violating transitions, flavor-blind quark or 
lepton observables, or at high energy colliders (see e.g. Ref. [iS]). 

2. For all spins, there is in the basis an operator involving only the SM Higgs fleld coupled to the dark state. 
Though speciflc models may not generate such operators, their presence would have two consequences. 
First, in general, the invisible state can no longer be naturally massless, though it can be very light, 
since a mass shift arises after the electroweak symmetry breaking. Second, the NP scale should always 
be greater than the electroweak scale, A > w « 246 GeV, otherwise these corrections grow unchecked 
as powers of H^H are inserted (this is actually true for all our effective operators). It should be noted 
though that for spin 1 and spin 3/2 particles, enforcing a dark gauge invariance on the couplings to SM 
fields explicitly forbids such mass terms. 

3. The contribution to FCNC transitions of the flavor-blind operators, i.e. either those involving SM gauge 
and/or Higgs flelds or those involving SM fermion flelds of the same flavor, has been clarified. Specifically, 
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for scales A much higher than the electroweak scale, we pointed out that it is always advantageous to dress 
flavor-blind operators at the low scale with a W exchange (see Table [5]) . Indeed, such a low-scale GIM 
breaking does not necessitate additional Higgs tadpoles, which would each bring in a 1/A suppression. 

4. For spin 1 and 3/2 dark particles, special care was devoted to maintaining a sensible massless limit, or 
at least to interpret the seemingly divergent limit. Indeed, the flavor-changing neutral currents are not 
conserved in general, so that the l/my (l/rn|,) term of the polarization (spin) sum is not projected out 
in physical observables. This is particularly relevant for massive vector states, for which renormalizable 
couplings to non-conserved quark currents can be constructed. Several mechanisms were discussed, and 
the corresponding NP scales derived from experimental bounds on the rare decay branching ratios were 
compared. As shown in Figs. S] and O these scales strongly depend on the assumed dark sector dynamics. 

5. All the leading operators producing a single dark fermion, whether of spin 1/2 or 3/2, violate either 
baryon (B) or lepton (C) number, but not both simultaneously, and have dimensions smaller or equal 
to that of the FCNC operators. By contrast, most of the AB and AC violating operators involving a 
dark vector or scalar particle directly derive from the dimension-six AB = AC = 1 Weinberg operators, 
or from the dimension-five AC = 2 operator. The former are negligible when they induce proton decay, 
i.e. when m^y < rrip, but could induce exotic B decays into an odd number of baryons plus missing 
energy when niB > > nip, since then proton decay is kinematically forbidden. The latter AC = 2 
operators do not induce the quark flavor transitions, hence have a negligible impact on rare decays. 

6. To implement the SU{2)l <S) U{1)y invariance, the FCNC operators are constructed in terms of the SM 
chiral fermions. Though theoretically sound, and particularly convenient to implement the MFV flavor 
restrictions, such a basis is not convenient phenomenologically because many of these operators interfere in 
physical observables. So, our bounds are always derived using the alternative basis of operators obtained 
by projecting on currents of definite C and P. This minimizes interference terms, since in most cases 
these currents produce different final states. 

7. Although many examples of NP models involving new light states were mentioned, no attempt was made 
to precisely match them onto our basis of operators. Indeed, this would require dwelling into the detailed 
dynamics and parameters of each model, and would bring us too far from our main objectives, which 
were to construct the most general basis and constrain its operators from rare FCNC decays. 

Concerning the rare FCNC decays of the K and B mesons, let us remind that provided the non-standard 
light states are neutral and sufficiently long-lived, they would show up as missing energy. Those modes cannot 
be experimentally distinguished from the SM processes producing a neutrino pair in the final state. However, 
these SM processes are very suppressed, and thus in principle, the rare decays could permit to identify even 
tiny NP effects. The main outcomes of our detailed phenomenological analysis of such effects are: 

8. First, we stressed the importance of including the correct kincmatical dependences for probing NP op- 
erators. This is crucial because experimentally, the rare decay modes with missing energy do not allow 
for a complete kincmatical reconstruction, and require aggressive background suppressions. In practice, 
most experimental analyses implicitly assume at various stages that the differential rates have the shapes 
predicted by the SM, and this seeps through down to the final bounds on the branching ratios. Note, 
importantly, that these dependences are not always accounted for by simply enforcing the various exper- 
imental kincmatical cuts. For these reasons, the specific kincmatical dependences of each NP effect may 
have to be implemented by the experimentalists (those are detailed in App. A for K decays, and App. B 
for B decays). This provision has to be kept in mind when interpreting our bounds. 

9. In the K sector, the sensitivities of the Kl X , K ttX, Kl — 7A', and K ttttX channels 
were compared, see Tables [3] and El with X a single or a pair of dark particles. The two-body K X 
mode turns out to be the most sensitive, though for a very limited number of operators, but it is also 
the most difficult to deal with experimentally. At the other extreme, the K — ttttX channels are 
sensitive to nearly all possible operators, but do not appear competitive given their phase-space and chiral 
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suppressions. This leaves the K — )• nX and Kl — >■ channels, whose sensitivities to NP operators are 
in general comparable. Still, it should be noted that the latter, not yet considered experimentally, has 
some advantages. First, its SM contribution Kl ^vv is at the 10"^'^ level (see App. lA.S.ip . and thus 
cannot obscure even a tiny NP contribution. Second, for massive flavor-blind dark vector bosons or for 
millicharged fermions, as derived e.g. from a dark U{1) kinetically mixed with U{1)y [28], the — ^X 
mode is significantly superior to if — > vrX, whose relevant matrix elements vanish at leading order in 
Chiral Perturbation Theory. 

10. In the B sector, the sensitivities of the Bg^ —t' X , B {K, K*)X , and B — {tt,p)X decay channels 
were compared, see Tables H] and El As for K, the fully invisible decays are both the most sensitive and 
the most difficult to probe experimentally. The main new feature compared to the K sector, besides an 
extended kinematical range allowed for uix, is that the modes with two light mesons in the final states 
can resonate, so that B — )> K*X, and B — > pX are competitive. This is particularly interesting since 
these modes are sensitive to all the quark currents but the scalar bs and bd. It should be said also that 
the present sensitivity of 6 — ^ s and b ^ d decays, in terms of NP scales, is very similar. So, if the NP 
operators do not follow an MFV-like scaling, a small NP effect could be easier to identify in the latter. 

11. The relative sensitivity of K and B decays was also compared. As expected, if the flavor structures of 
the NP operators involving X are generic, K decays are far more sensitive than B decays. However, it is 
well-known that in the visible sector, generic NP flavor structures are at odds with current experimental 
constraints. If MFV is imposed on both the visible and dark sector operators, the constraints from B 
decays become often tighter than from K decays (see Fig. [3]), especially for chirality flipping currents 
5^(1, 75, (t'^'')^"', relatively suppressed by ms/mi,, and for low-dimensional operators. Indeed, the impact 
of MFV on the scale A for an operator of dimension n decreases as n increases, since it is approximatively 
given by Amfv/A w (VjjVtj)^/^""''^ for the d'^ d^ transitions, and with the CKM coefficients given in 
Eq. ([5]). Note, however, that n cannot be too large, since rare decay constraints give A < v when n > 8, 
see Table ([I]). In other words, for A > w, the impact of such operators on the rare decays is beyond reach. 

12. The AB and AC operators have low dimensions only for X ~ tp or "if. The AB = 1 operators can only 
be probed with specific searches in B decays involving an odd number of baryons plus missing energy 
in the final state, and should certainly be included in future experimental programs. For the AC = 1 
effects, the low-dimensional operators are accessible only for X = i/j, which contribute to P+ £^tp 
(P = K,D,B), and would thus apparently enhance the purely leptonic transitions l^v. If the 
flavor structure of these NP operators is non-universal, this could resolve the persistent discrepancy in 
B ^ TV while remaining consistent with the B — > (e, p,)v bounds, as well as, if < rriK, with the tight 
lepton universality constraint derived from Kn decays, see Eq. pT|l . 

In conclusion, the presence of a light invisible state weakly coupled to SM particles is not only far from 
excluded, but is even compelling in many NP models. To flnd such states, a host of experimental facilities 
currently available or in planning are called in, from high-intensity meson and lepton factories to high energy 
colliders, neutrino detectors, earth or space-based direct or indirect dark matter searches, high intensity lasers,... 
In this big picture, the very rare FCNC decays of the K and B mesons, with their unique sensitivities and 
kinematical ranges, could play a crucial role in the very near future thanks to the leap in luminosity expected 
at the next generation of experiments, namely NA62 at CERN and KOTO at J-Parc dedicated to these K 
decays, and Super-B in Italy and Belle II at KEK aiming for the B decays. 
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A Differential rates for K decays 

A.l Experimental observables in rare K decays 

Let us start by reviewing the kinematics and current experimental limits for the various K decays induced by 
neutral currents and involving missing energy. 



A. 1.1 K 7r+ missing energy 

When the missing energy consists of two invisible particles, the differential rate depends only on the invariant 
mass of these particles, z = q'^/m^j^, or equivalently, on the pion momentum Pj^ = \v-n\l"n^K = \/^/2, with 
A = A(l,z,r^) defined in Eq. r^r = mT^/mK- The phase-space integral is then 



dz^[z]^ ,7 (66) 



with rx — nix/mK- 

In the SM, the only available invisible particles are the neutrinos. The SM spectrum for A'+ — > n^iyi? 
and Kl — > tt'^vv then derives entirely from the vector current matrix element {Tr\sj^d\K), and involves the 
corresponding form- factor (see Eq. (I5T|) in the next Section) slopes A'^ and A" : 



,r-(l-r,)^ 



/f - (0) 



^'""^ dzA3/2fl + AV-^ + A^|_) , (67) 



where i = +,0. Translated in terms of the pion momentum, i.e. using z{P-^) = 1 + — 2 i/r^ + P^ , this 
becomes 



•/ri + p; V ^ ^ 2'-; 

The slopes A^ and A" are conventionally normalized by the charged pion mass, and are equal for the and 
decays to an excellent approximation (see Ref. 59 ). They are extracted from K13 decays as 

A^ = rx (24.82 ± 1.10) • 10"^ , X'l = rx (1.64 ± 0.44) • 10"^ , rx = 0.990(5) , (69) 

with rx a rescaling factor accounting for the A'*+ — K*'^ mass difference. These (highly correlated) errors are 
negligible compared to the experimental and theoretical errors on the integrated rate, and are neglected in 
Table H 

Currently, only the charged decay has been observed at Brookhaven [60j . in two momentum windows 
separated by the — > tt+tt'' peak, and with the lower end corresponding to the K — ?> tttttt threshold (see 
Tabic [5] and Fig. 15]). The proposed charged K experiment at J-Parc would use the region above the — >■ tt+tt^ 
peak [61,, while the two windows planned at NA62 [62 and proposed at Fermilab [SJ ^"^^ similar. Note that 
these experiments use very different techniques (stopped vs. in flight), but in both cases, the momentum of 
the initial and final charged particles are in principle accessible. It is important to stress that not only the 
combination of the measurements done for each specific window (see Table ^ assumes the SM spectrum, but 
also that within each window. 

For the neutral mode, the Ks — n^viy mode is CP-conserving but difficult to access given the very short 
Ks lifetime, so we concentrates on the CP-violating Kl 'n^vv mode. The best limit [B3] 

B{Kl -> -K^vv) < 2.6 • 10^^ , (70) 

was obtained by the E391a experiment at KEK, and will be further improved using the same techniques at 
J-Parc |61| . In these experiments, the momentum is not fixed. So, the high hermeticity of the detector 
is essential to ensure sufficient suppression of the backgrounds. Though less effective in this case, kinematical 
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\p.\ {MeV) 


z 




SM 


Extrapolated total 


[211 , 229] 
[ 140 , 195] 


[ 0.000 , 0.062 ] 
[0.116 , 0.289] 


27.6% 
39.7% 


0.228(18) 
0.328(25) 


1 47+i-3U 

-0.89 

7.89^^;?^ 


Combined (Total BR) 






0.825(64) 


1 70+1.15 

^- '«J-1.05 



Table 9: Experimental measurement of — > tt+j/j/ and SM prediction within each momentum window |60j . 
in units of 10^^°. 




P;il ^'^^^^ 

Figure 6: The experimental windows in 7r+ momentum used for controlling backgrounds in the — >■ tt'^vD 
measurements, with the seven events seen at Brookhaven. The SM spectrum corresponds to a vector coupling 
sjfj,d X PlJ^vl, and is implicitly implied in computing the branching ratios from the events. 

cuts are still useful. In particular, the transverse momentum Pt of the reconstructed tt" is required to be large, 
between 120 and 240 MeV. This does not cut away the background from Kl tt^tt^, but rather ensures that 
the two extra photons have high momentum, and are thus difficult to miss. Since the momentum spectrum of 
the 7r° cannot be directly measured at KEK or J-Parc, and since the SM decay spectrum is implicitly assumed 
in the analysis, it is far from immediate to translate the current limit (|70p into bounds on non-standard currents 
involving other types of invisible particles. 

So, for both the charged and neutral modes, it is not currently possible to deconvolute the SM spectrum 
from the experimental numbers. To proceed and derive the bounds quoted in the text, we require that the 
predicted branching ratio for the production of new invisible states does not exceed 10~^° when integrated 
over the momentum windows of the charged mode. This is a rather loose approach, which could significantly 
underestimate the experimental reach in case the spectrum is very different than the SM one. To illustrate 
this, note that the bounds for two-body decays are already slightly tighter [65] . 

B{K+ TT+X°) < 0.73 • 10^1" {mx = 0) . (71) 

To improve our naive bounds, either the specific modulations of the spectrum in the presence of NP have to be 
included throughout the experimental analysifH, or the true momentum spectrum must be measured (maybe 
using TOE techniques for the neutral mode [53]). Einally, independently of the NP spectrum, it should be noted 
that the sensitivity to the K — )■ ttX{X) processes is ultimately bounded at around 10""'^^ by the theoretical 
error on the SM predictions for the K — > ttvP branching ratios. 

^Brookhaven extracted bounds for purely vector (SM), scalar, or tensor currents, but not for a combination of the SM plus 
non-standard interactions. Since we cannot turn off the SM rate, this is not directly useful for our purpose. 
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A. 1.2 K — > 7171+ missing energy 

The phase-space integration for the K — > 7i{Ki)7t{K2)X [q] decays is 



dy i- , (72) 

4^2 ay 



with y — [Ki + K2)^/mj^ the invariant mass of the pion pair. Compared to K ^ 7tX, these modes are 
suppressed by the smaller hadronic matrix elements and by phase-space. Further, the kinematical range is 
much reduced. Currently, the best limits are (see the respective papers for different mx values) 

BiK+ ^Tt+TT^X") KA-lQ-^mx ^50MeV] gS], (73a) 

B{Kl 7r°7r"X") < 7 • 10"^ [mx = 50 MeV] [57], (73b) 

and are thus very far from Eq. ([7T|). Note that the hadronic matrix elements (tt+tt" 

and {TT~^Tr~\sTd\K^) are related in the isospin limit, see Eq. (IM]) below, so that these experimental constraints 
suffice to completely bound the K — ^ ttttX system. 

The K {P) — 7> 7t{Ki)'k{K2)X(pi)X{p2) decays are similarly suppressed, and the experimental information 
is less precise. The phase-space integrals reduce to that over the invariant mass of the invisible pair = 
{Pi +^2)^ = zm^ and of the pion pair = {Ki + K2Y = ymj^, over the range 

dz dy—-. (74) 

r2 J4,.2 dydz 

Here again, the SM spectrum critically enters and the current experimental bound 

B{K+ TT+TT^iyD) < 4.3 • 10"^ [66, , (75a) 

B{Kl 7T°7T°iyD) < 8.7 • 10"'' [67] , (75b) 

cannot immediately be translated into bounds for invisible particles of a different type. 

Though probably optimistic given the limited phase-space and complicated signatures, we assume bounds 
of 10^^° on each of these K — >■ ttttX and K — ?> ttttXX modes are achievable to derive the numbers in Tables [3] 
and [5] In any case, if the bounds are different, it is a simple matter to rescale the numbers accordingly. Further, 
using the same 10""'^° branching ratio bounds as for the K 7tX{X) modes permits to clearly illustrate the 
reduced sensitivity of the K — )■ 7r7rX(X) channels. 



A. 1.3 K 7+ missing energy 

Compared to K ^ ttttX and K ttttXX, the modes K (P) -> 'y{k)X{pi)X{p2) and K (P) ^ j{k)X(T) are 
less suppressed and could offer simpler experimental signatures. The phase-space integral for the three-body 
decay is 

^ixx = dz — , (76) 

with z = (pi -\-p2Y lm?j^ the invariant mass of the invisible pair. Despite their theoretical sensitivity, there is 
currently no experimental limit on these modes. So, for now, we assume that the next generation of experiments 
will reach B{Kl — > ^X) < 10^^°. Note that with about 10^^ — 10^'^ Kl decays, as required to measure 
Kl — > 7r^vi', this may be pessimistic. Further, while the K — > 7rh'iy processes ultimately limits the sensitivity 
to K ^ 7rX{X) at a few 10"^^ given the current theoretical errors, the SM rate for Kl — > Ji^i' is at the 10"^'^ 
level, so bounds at or even below that level are in principle achievable. 

Other modes with photon and missing energy will not be considered, as the K — > titt + J + X processes are 
either too suppressed and difficult to access experimentally when fully neutral (given the many photons from 
the tt'^s), or superseded by the non-radiative processes K — titt + X when some mesons are charged (since 
the photon ot K mr + j + X is essentially a bremsstrahlung radiation off the charged meson [68l|69] , the 
amplitude for if — >• n7r -I- 7 -|- X is actually driven by that ioi K ^ mr + X [70]'). 
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A. 1.4 K — > missing energy 

The simplest decays are those where the or Kg simply disappear. Though difBcult to probe experimentally, 
the simpler matrix element together with the minimal number of final state particles strongly enhance their 
sensitivity to NP effects. In the case of the 7r° — > XX process, the best bound is [71] 

e(7r° ^ XX) < 0.27 • 10"^ . (77) 

But this measurement is actually a by-product of the study of the — > ti^XX decay, since a bound on 
A'+ — > Ti^XX indirectly constrains A'+ — > 7r+7r"[— > XX]. Doing the same for Kl — > XX would require very 
tight bounds on some B oi D decays with missing energy, well beyond current capabilities (see Sec. IB.1|) . 
Alternatively, a direct bound on Kl XX may be obtained from (p factories, where the other K can be 
tagged. In deriving NP scales in the text, we will use B{Kl XX) < 10^^" to simplify the comparison with 
the other modes, but it should be kept in mind that such a bound appears extremely challenging. 



A. 2 Matrix elements for K decays 

In the K sector, the quark currents are represented within Chiral Perturbation Theory (ChPT) 52 . For 
simplicity, only the leading chiral order is kept. Specifically, the vector and axial- vector currents start at 0{p): 



l^^.J (^Df^U^Uy'', qI^^^^qi^^—{D^^UU^)", (78) 



2 V 2 

with, at leading order, F = F^^ = 92.4 MeV. Thanks to the QED gauge invariance, there is no unknown 
low-energy constant in these currents. The scalar and pseudoscalar currents start at 0{p^): 

^Wr^-^B^U", ql^qi^-^B^U^" , (79) 
with the low-energy constant Bq related to the quark masses, 

B.^^^ (1-^^=4.6(8), (80) 
■niK rus 

with r^r = mTr/niK, and using ms(2 GeV) = 100 ± 20 MeV (so when deriving bounds on q/A", Ci = Ci{2 GeV) 
is understood) [H]. The scalar and vector currents are related by the EOM. Specifically, the most general 
matrix elements for the scalar or vector K ir transitions have the form (z = (f' /fn^K-i Q = P ^ ^) 

{n [K) \sd\K (P)) ~ "^^"""' /o {z) , {7T {K) Is^dlK (P)) ^ (P + i^)^ /+ (z) + (P - K)^ /_ (z) , (81) 
nis - rud 

up to some simple Clebsch-Gordan coefficient. Taking the divergence of the vector current produces q^f- (z) — 
{m\ — m^) (/o (z) — /+ (z)). So, at the leading chiral order, /+^o (z) = 1 and /_ (z) — 0. Refinements are 
only needed for a precise prediction of the SM rates, but are not numerically relevant for the bounds on the 
production of new invisible states. Note that in practice, the scalar and vector currents do not need to be 
parametrized as external couplings, but can be directly introduced through the ChPT source terms. Doing this 
using the leading 0{p'^) Lagrangian reproduces Eq. (175)) and Eq. ([7^ . We do not consider the next-to- leading 
O(p^) meson loops and local terms, except for the odd-parity contact interactions obtained by introducing 
the vector and axial vector sources in the 0{p'^) anomalous Wess-Zumino-Witten (WZW) action [52]. Indeed, 
owing to their opposite parity, these interactions drive the leading order contributions for some amplitudes. 

Finally, from Lorentz and chiral symmetry, combined with parity and charge conjugation (valid for the 
strong interactions), the most general parametrization for the tensor current starts at 0{p^), where it is given 
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by m 

+ ^a'^iiFj:, - ^Fl:^)U^ + U\F^^ - ^F^^))" , (82a) 

+ ^a'^iU{F^, + zF^^J + (F^f, + iF;i)Uy" . (82b) 

Two new low-energy constants qt and aip occur, for which we use the Lattice estimates (see the discussion in 
Ref. i5U) 

Bt{2 GeV) ^ 2mKaT = 1.21(12) [72^ , 5^(2 GeV) = 2F^aT = 0.6(2) [73] . (83) 

Note that a more recent lattice estimate Bt — 0.65(2) [71] is two times smaller, and thus suppresses the 
sensitivity of K decays to the tensor currents. Still, in terms of the NP scales A given in the Tables [3] and [S] 
the precise value of Bt is not that relevant at present since these numbers are to be understood as order of 
magnitude estimates. 

So, altogether, and defining F = + cpj^ + cyj^ + Cyi7^75 + cto^i, + CfUp^^^, the matrix elements in the 
isospin limit and to the leading chiral order are|^ 

V2(0|srd|X°(P)) = 2iF^ {-Bocp + PPcA) , (84a) 

V2(7(fc,a)|srd|7^0(P)> = -jeF„{a'j.Cf {k^^g'^'' - k^^g^P) - la'rCTe^P'^Pkp + ^^e^PP^kpP,) , (84b) 

V2{Tr" {K)\sTd\K"iP)) = -{7r+ {K)\sTd\K+ (P)) 

= -Bocs + (P^ + KP)cv + 2aT (ictK^^'P"^ - CfE^^P^KpF,^ , (84c) 

and 

-(7r+(/vi)7r"(i^2)|srd|/^+(F)) =M-{K^ vrvr) , V2{TT°{Ki)TT°{K2)\srd\K'^{P)) = M+ {K tttt) , 

V2{7r+ {Ki)'K- {K2)\srd\K° (P)) = M+{K ^ 777:) +M-{K ^ t:?:) , (85a) 

with 

M,iK ^ nn) ^ - f (^^ + ^^f^^) + f {^f^l'T^' ~ .c.e--T_,AV.) , (86a) 

(86b) 
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where K± = K2 ± Ki, T± = P ± and X^^V^ = X^'Y'' - X''Y^'. Terms proportional to the number 
of QCD colors, Nc — 3, come from the WZW action. The m\ — denominators arise from the kaon pole 
topologies, K — ^ ttttK^ followed by (from Eq. (IMl)). Note that M + {K ^ tttt) is even under Ki ^ K2, 

while A4^(K — tttt) is odd, hence these amplitudes describe two-pion states with even and odd orbital angular 
momentum, respectively. The tt^tt" state is purely even due to Bose statistics, while the tt+tt^ state has total 
isospin one, hence is purely odd. 

Only the tensor currents contribute to the Ai{K 7) amplitude at tree- level. So, it may seem that 
together with the O(p^) WZW amplitude, we should include also the even-parity meson loops along with their 



®A number of sign conventions are implicitly defined by these equations. We closely follow the ChPT conventions of |51II52| 
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counterterms. However, for neutral current sources, i.e., in terms of Gell-Mann matrices, for v^^a^, s,p ~ 
Ae ± iAy, the only allowed even-parity X — !■ 7 matrix elements vanish when the photon is on-shell 



72^/27^2^V 



q^T" -q^q-T 



$(g^m^)-^((7^m^) 
g2 



where T = P — q, and $(9^, m^) the loop functions occurring for K — ttj* (see e.g. Ref. [55]), defined in terms 
of the standard scalar one-loop integrals as 

$(g^m2) = 3(g2 -4m2)Bo(g^m^m2) + Um^BoiO^m^^m^) - 2q^ . (88) 

Contrary to JsT — )• tt'j* , the FCNC matrix elements are finite since the UV divergences cancel in the difference 
between the and tt^ loop contributions, and there are no counterterms. Though in principle, the K — > 
£'^£~X modes could thus offer precise probes, their rates are far too suppressed by a, the loop factors, and the 
cancellation between the and tt^ loops. So, only K — >■ ■fX will be considered here, which is thus induced 
exclusively by tensor currents and 0{p*) anomalous interactions. 

Finally, the operator basis also includes vector and scalar currents with a covariant derivative. Though 
these operators are never retained in deriving bounds on the new physics scale in the main text, for the sake of 
completeness, let us nevertheless write down the relevant chiral realizations. For the vector currents, extending 
the analysis of Ref. [75^ , we write 

q'^Da^j^PLq' = i^id^U^d^U + d^d^U^U - ^,gap(xtc/ + U^x) + avie^^p,df'U^d''Uy" , (89a) 

q'^Da^j.^PRq'' = i^{d^Udo,U^ + do^d^UU^ - \gc.^.{}(U^ + Ux^) - avie^^^^d^Ud^'U^)" , (89b) 

q'^pL^aq' = -i^idcU^d^U + U^d^d^U ~ \9a.^{U^x + X^U) + ayie^^fs^d^U^d'^Uy' , (89c) 
4 o 

q'l^PR^cq' = -i^idcUd^U^ + Ud^d^U^ - Ig^^iUx^ + xC/t) - aviSc^^p^d^Ud^U^y" . (89d) 

4 8 

Most of the terms are fixed by taking divergences and imposing the EOM, but for the constant ay, a priori of 
0(1). For the scalar currents, the chiral realizations start at 0{p), 

q^PR^^q' - q'Pi.'D.^q'' - -^Bo{D^Uy" , (90a) 

q'PLD^q-' - q'PLD,,q-' = ^^Bo{D^U^y' . (90b) 
These currents are completely fixed by imposing parity and charge conjugation, together with 

fPL,R^aq' + f^aPLMq' - d^{q-^PLMq') ■ (91) 

Note, however, that the above chiral representations lead to (f (J^ a — ^ a)PL,B.q'' — instead of {inj — 

m'j)q^ Pl^bM'^ ■ This is because while q^(ia + ^ a)PLMff' is of 0{p), the difference q^(l^a — ^ a)PLM.q^ is 
actually of 0{p'^) since quark masses are 0{p^). So, terms at that order would be required to get a correct 
divergence. 

A. 2.1 Decay rates in the isospin limit 

The strong matrix elements (IMH55|) are derived in the isospin limit. As a result, all the differential decay rates 
can be reconstructed entirely from those of the Kl ~ K2 ■ The contributions coming from the eKi piece of the 
Kl, suppressed by £ ~ 2 • 1Q~^, are neglected here. 
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For the K — (7)^ modes, the Ks ~ Ki rates are obtained from those for Kl ~ K2 by interchanging the 
real and imaginary parts of the couphngs x = fi, gi, hf. 

dV {Ks ^ (7)^) dT (Kl ^ (7)^) 



dz 



dz 



o 5R(x)] 



(92) 



The K^ — >■ TT+X decay rates are proportional to the sum r{Ks — > tt'^X) + T{Kl ^ tt^X), hence are obtained 
from T{Kl — tt'^X) through the substitutions {x, y = fi, gi, hi) 



dr {Ks Tr"X) _ dT {Kl tt°X) 



dz dz 
dr{K+^Tr+X) _ dT{KL^TT°X) r 3(a;)2,5ft(x)2 ^ Ixp , 



dz 



dz 



^{x)^{y),^{x)^{y) 5R(xy* 



(93a) 
(93b) 



Finally, the whole set of K ttttX decay rates can be reconstructed from the Kl 'n'^'K~X rate as follows. 
Denoting T{Kl — >■ tt^tt^X)^ ^ \M±{Kl — 7r7r)p from Eqs. ([Ml [55]) . and noting that these two amplitudes 
do not interfere, we get 

d^T {Kl -K^-K^X) 1 rf^r [Kl ^ tt+tt-^)^ 
2 dzdy ' 



dzdy 

d^T {Kl ^ n+n-X) _ d^T {Kl ^ ir+n-X)^ d^T {Kl ^ tt+tt-X) 



dzdy 



dzdy 



d^T {K+ 7r+7r"X) _ d^T {Kl ^ TT+7r-X)_ 



dzdy 



dzdy 



dzdy 

^{xf,^{x)^ -> \x\^ , 
^{x)^{y),^{xmy)^^{xy*) 



d^T {Ks ^ 7r"7r"X) _ d^T {Kl ^ TT°n'>X) 



dzdy 



dzdy 



(iiTs -> TT+TT-X) _ (i^i -> TT+TT-X) 



;3(a;) o 3fi(a;)] 



(94a) 
(94b) 
(94c) 
(94d) 
(94e) 



dzdy dzdy 

In the following sections, the differential rates are given for the various scenarios adopting the notations of 
Eq. dMl) for K ttXX modes, Eq. ^ for K ttttXX modes, Eq. ([72]) for K ttttX modes, and finally, 
Eq. (I75|) for K — > -fXX modes. For K — )• ttX and -fC — jX, the total integrated rate is directly written down. 
Given their regular occurrence, let us also introduce specific notations for the usual kinematical functions. 
First, 

= A(l,a, /3) , X{a,b,c) = a'^ + b'^ + ~ 2{ab + ac + be) , (95) 

with a,/3 standing for the reduced variable y, z, or the reduced mass r^r, rx (in which case we simply denote 
— TT,X). Similarly, we define 



Pl = l- 4rl/y , /3f - 1 - 44/z , /3: 



1 - 4rf 



(96) 



A. 3 Spin 1/2 invisible particles in the final states 

The rate for the fully invisible decay is: 



4. 



TO] 



r {Kl ^ = {Ii3(/ps)2 +Iig?(/pp)2 +I25R(/aa)' +Xi25R(/aa)5R(/pp)} 



Fn = Xi = I( = , I2 = 4rL Z]2 = -4r,„B 



2^TOK ' -1 - - ^' - ■ 

The differential rates for the decays into a pion plus invisibles are 

^ {Kl ^ = ^FJ {^i3(/tt)' + :^i'5i(/fT)' + J2^{Ivv)^ + Jl2^{fTT)^{fvv) 

+J2HfvA)^ + Jsmss)^ + JMfsp)^ + J23^{fvA)^{fsp)} , 



(97) 



(98) 
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with the normalization and kinematical functions 

J23 = -ir^Bs (1 - tI) , ^4 = J5 = ^B^s , Btz = ^B^Xz. ■ (99) 

In the massless limit, this expression simplifies a lot because the fermion helicity states do not mix. The 
interference terms drop out while the parity of the current becomes indistinguishable; J7i — J7^', i ~ 1,2,3. 
The decays K — >■ mr^^p also receive contributions from all the currents: 

+-^75? {fvvf + ^^3? ifvAf + -^275? (/tt) 5R (/yv)} , (100) 
'^'^ (i^i ^ vr+Tr-V^V) , = ^rj- {-F45(/P5)' + -Fl3fi(/pp)' + HmjAA)^ + ^ismAAMfpp) 



dzdy ^ " ^^^+ A4 

+-F53?(/Ay)' + Te^ifTTf + J^e'^iffrf + J^BemtrMfAv)} , (101) 

with the normalization and kinematical functions 

= ?||^, Ar = /3^(y(l - yf + ^(4y + .)), A. = /3^(,z + ^), 
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2(1^' " 2 (l^z)2 ' = ^6 = 2//?,'A,., J-i = -fiXiPi + i84/z)Xe), 



?2 



= ^(Ai(l + &rl/z) - {8rl/z)Xe) , ^"2,5 = A2^5(l + 2rl/z), T'^^^ = T'l^ + l^Ag^, -^2,5 = As.s/?: 



2 



— V^viV-L ^ i/)/ ^1 V"' ^/ , -I 2,5 — '>2,5V-^ ^' V/ 2,5 — 2.5 ^ 2,5 ~ ^'2,5/^2 , 

^ 77" TT" _ "'^ \ d2 tt _ d2 ^Ayz „9 „9 ZA„ 

I „\ 77' — 77"fl2 7-'/ 



-^3,4 - -7^3,4/^2> -^3,4 - 9-^3,4' > -^M - 9A3^4B3 , - B^—^/Sj , J"g - B^—^{1 + 8r|/z) 

J-7 = J-^' (1 + 2rl/z) , = , j;' = ^Ag^^^M/, -^27 = -r^BrXeAwzw , 

Tx2 = ^BtPI {Xy, + 12y(l - y)) , ^"23.45 = -r^A3,4S„ J-ge = -r^SrA^, , Awzw = ^T^rf • (102) 

Again, the interference terms drop out and T[ — i = 1, 6, when — > 0. 

The rate Kl — >■ ^ipip is driven either by the anomalous vertices at 0{p'^) for vector currents, or directly by 
the tensor currents, and has the differential rate 



° " 275^2 ' ^'1 = — (1 + 2r^/z), ^^2 = 3r^ByA^^ , 



c;3 = ^A^^, ^?4 = ^(1 + 2r^/z)4,, A,, = , (103) 

where A^c is the number of QCD colors. The O^p^) loops and normal-parity counterterms cancel out for the 
(axial-)vector and (pseudo-)scalar currents when the photon is on-shell, see Eq. (1571) . 
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A. 3.1 Standard Model rates 

The standard model rates are recovered by setting all the coefhcients to zero but for 
AGpa (Mz) y. 



'^LL 

A2 y/2 27r sm^ 9w 



(104) 



Numerically, Xt = 1.465(16) [Z^, Pc = 0.372(15) [77], SPu,c = 0.04(2) [7S] (with A = 0.2255), so that for each 

^ — C, /i, T, 

y^ = 27rsin2 6'vi/ x [4.84(22) - il.359(96)] x lO"'' . 
The full set of differential rates is, in the isospin limit, 



(105) 



TiKL 







X / dzXli^ X 3(y,)2 



TT n vv] 



dzdy 



P.X 



•'K 



647r2F2 



X ^{y,f 



T{K, 



n^TT vv 



) _ = Ti^p X y dzdy 



r {Kl ^vv) = Ti^p X / dz 



192tt^F3 



arnj^zll — z)^ 



^{y,f X {X + 12yz) + n{y,f x ^ 



yzXyzm%N, 



t21 



727r4Fi 



X 3?(y^ 



(106a) 
(106b) 

(106c) 

(106d) 

(106e) 



where Fi^p ~ G pa (Mz)'^ m% / {256tt^ sin^ 9w), and the ranges for the phase-space integrals are given in 
Eqs. ((66l [74l [76l) . Note that -RTi, -> 71^/? is purely CP-conserving because the parity even matrix element 
vanishes at 0{p'^), see Eq. ([57)1 . 

The corresponding branching ratios are 



B {Kl -> ttVp) = 2.3 x 10~", B {K+ ir+vv) = 7.6 • 10"", 
B (Kl ^ TT+TT-vv) = 1 • 10"", B {Kl ^ TrV^i^P) = 6 • 10"^"^, B {K+ tt+tt^vv) = 5 • 10"" 
3.4 X 10"", B (Kl ^ e+e-vv) = 2 ■ 10"", B (Kl fi+fJ.-v9) = 7 ■ 10"" 



B{Kl 



•yvv) 



(107a) 
(107b) 
(107c) 



Adopting the usual chiral counting, the typical error on these LO estimates is expected to be of about 30% at the 
amplitude level. Note, in this respect, that the K — > ttvv rates given above are just indicative, as higher order 
corrections as well as isospin-breaking effects are known and more precise estimates have been obtained [59] . 
For K — >■ TTTTz/P, the anomalous term gives negligible percent-level contributions, so that '^(yu) < ^(yv) 
implies B {Kl tt+tt^vv) « 2B [Kl tt'^tt^vD) > B [K+ tt+tt^vv), in fair agreement with Ref. [75]. For 
Kl jvv, previous estimates incorrectly rely on Ki2-y for the matrix elements [Hn]j hence included a parity- 
even contribution in contradiction with Eq. (1871) . Numerically, this is however without consequences since these 
contributions are CP-violating hence subleading for the rate. Finally, the rates for Kl — >■ t^l~vV are dominated 
by the Dalitz emission from the purely anomalous Kl — >■ ^vv. The rates from the even-parity contributions 
arising from the matrix elements Eq. (157)) are in the 10~^^ range. Note that we did not consider the tree-level 
process Kl — >■ W^+W^~[— i'^i'v^vi], which may actually be competitive given these strong suppressions. In 
any case, these rates are far too small to be accessible any time soon. 



A. 4 Spin invisible particles in the final states 

The rates for the production of a single invisible scalar from the H^{DQ)(p operator of Eq. ([2(1)) are 

dr 



F {Kl ^ vrV) - ^B^Mys? , 
F {Kl ^ 7</>) = , 



, {Kl TT+TT (/))_ 

— {Kl 7r+7r"(?!)) , 
ay ^ 



f J- ni^{gpf , 
f 5- H23(5p)' , 



(108) 
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with 

> 1/2 0,1/23 / \ 2 2 \ 

V mK Py\y rUj^ I y ^ Py\^y 

- ' = T024^' ^1 - - Y37| j = 3(l-r2)2^^ ' ^^^^^ 

The matrix element (715^1/^^,) occurring for — >■ 70 vanishes at 0{p'^). The rates corresponding to the 
derivative couphngs {s^^d)d^(j) and (57/^75^)9^0 are obtained through the replacement (1^^ . i.e. 

. rUsTrnd .,, l-r^TOK , , 1 tuk 
5s,p -> ~i9v,A 7 with ms-md = — -— and rn^ + ma = — 7— . (110) 

A _Ds A _Ds A 

Note that if simultaneously present, the gs,p and gv.A currents obviously interfere. 

The differential rates for two-scalar final states can be written in terms of the same kinematical functions 
as for fermionic final states, but for the obvious replacement — > everywhere. They are significantly 
simpler though, because the angular momentum of the two-scalar states is purely orbital. Specifically, 



^ {K, ^ .^.-U)^ = {-^^'^^(ffAv)^^ +^^5R(5..)^^ + ^^'5R(gP5)^^ 



T^TTTT ( 4 



[K^ ^ 700) = ?^?33(5vy)'^ , (111) 



with the kinematical quantities defined in Eqs. ([Wl IMl [1021 1103P . Note that if = 0, Bose statistics has to be 
enforced and these rates should be divided by two. 

A. 5 Spin 1 invisible particles in the final states 

From the single dark vector production from the operators of "H^^JI] and "H^qJII], Eqs. (l28t and (|34|. we find 
^ vroy : rm = Yl^^{eyf , tM = f Av.^3(/t)^ , 



V 



A2 



^ : ^ ^ f5-(H^3(e^)2 +H^$R(£v)2) , = fr-rf {H'/5R(/t)^ + H^'5(/f )'} , 



with the definitions in Eq. (|109p together with 



(112) 



f7 _ 2amK(l-4)^ ^/ _ _ ,.2 Vy + 12y7-v . V^l .2 



H'/ = f B^yAv,/?^ ' ^2 = {Pl^- - ¥ + 12/3,^2/(1 - y)^) , K - S^^Av^, . (114) 



Note that the T-L'^ term is finite when ry — ^ since it is induced by the WZW anomaly, while by construction, 
all the r["l are finite in that limit. For the two- vector modes induced by the scalar and pseudoscalar couplings 
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of n^fnat [11] ' the rates are 



r["i {Kl ^ vv) = 4ro^ {i.mhpp)' + ^{hps)') + 64i[^ihps)'} 
dm 



rf2p[II] _ 6 

(i^i ^ 7r+^-T/y)_ = 4rr^ {zJ-3(5(/ipp)2 + 5R(/JP5)') +34J-^'5R(/ip5)')} , 

{Kl n+TT-VV)^ = ^ {zMmPpf + ^(hpsf) + i4T'Mhpsf} , (115) 

and the kinematical quantities defined in Eqs. (I97l99lf 02lf 03^ but for — my everywhere. The matrix 
element {'ylsdlKi,) occurring for — > 'yVV vanishes at 0{p'^). Note that if the vector field is real, V = V, 
then these rates have to be divided by two to account for Bose statistics. 

A. 6 Spin 3/2 invisible particles in the final states 

Introducing the short_-hands = ^{fx), 3?x = A[ = (S/S^' - 6/3| +9)/18, and /3f = (9/3,^ _ 6/3f + 5)/18, 

the rate for i^L ^ vj/vj/ is; 

4 1 

X5 = Ii^^, Ip = l(/3^^, = ^2f^lp , IpA = Il2P'^ , (116) 

where Fq, Ii, Z2, and I12 are defined in Eq. (P7| . The differential rate for the decay into a pion plus 
invisibles is 

AT jjA r2 

— {Kl ^ ^"^'*) = ^rj^ { Js^TS + ^p^Jtp + Jt^It + Jg^^f s + ^f'^fr + ^p^I p (H^) 

— Jpg'^TP^fS + JxP^Tt'^tp + Jpf^TP^fT + Jst'^TS'^TT 

— J gf^ftf g^fj, + 3yv(JVT3'TT + JTvs^^TS + iJypQ'fp) 

+Jss5i|s + Jsp3|p + Jy3yy + Ja'^Ia + ^Ap3y^3sp} , (118) 
with the normalization and kinematical functions 



16 t'^ I 10?'^ 



Jp — '^*2"g^ I ^ + \ ^■^PS — '^l~g^"^^-P — ^■^P^'^pf — g^l I 1 ^ I T^Sf — "^l"g^ I ^ 

JsT = 4Js, Jss = Js/^" , Jsp — JiP'z^ Jv = J2P'z — p^A^7r(3 — /3 ),Ja — Ji/S'z — T^—J^i, Jap = Ji2P'z, 



JVS = ^Jl2l3^,Jvp ^ ^Jl2 ( 1 + ~ + ~ ) ' '^^i^ = 9'^i2 I 1 - + I , (119) 



in addition to those in Eq. (|99p . The mode with a single photon with missing energy is 

rir 72 

- {Kl -> 7**) - AT r^i^X {^^^(^T5 + 4s) + Spi^'rP + ^^Ip) + eT(5R|T + 4t) 

+ fop(3?Tp3?f 5 + '^fs'^Tp) + QstO^TP^TT + '^fp'^frp) + Gpri^fp^TT + ^TP^f t) 

+ayKvy + Ga^Ia - '^vviGvT'^TT + Gvp'^fp + ^ys3?Ts)} , (120) 
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with the definitions in Eq. (|103p together with 



Finally, with the many possible interferences among the tensor currents and the complicated kinematical 
functions, the differential rates for Ki^ — >• tttt^'^' are too cumbersome to be given here. 



B Differential rates for B decays 

B.l Experimental observables in rare B decays 

Let us start by reviewing the kinematics and the current experimental limits for the various B decays induced 
by neutral currents and involving missing energy. For many of these modes, the kinematics and phase-space 
integrals are similar to that of the corresponding K decays, so we refer to Appendix lA.ll for the explicit 
expressions. 

Specifically, the observable three-body differential distributions for the B — > HX modes, H = tt, p, K, K* 
and X = -ipip, 4){4>), V(V), ^^E*, can be written in terms of the (reduced) invariant mass z = /m% of the invisible 
particles, or equivalently in terms of the H momentum in the B rest-frame (|p_ff | = 7715^^/2(1, r^)/2) or the 
missing energy {f] = mB{l + z — rjf)/2), as explicitly written down in Eq. for the K nX case. 

B.1.1 b ^ s+ missing energy 

For these transitions, we can obtain bounds on the production of new invisible states by using the existing 
experimental data from Babar [81][82] and Belle [83] searches for B — >■ K^*^i^P decays. In particular, the 
presently most stringent bounds are 

B{B+ K*+vv) < 8 • 10"^ [81] , (122) 
B{B+ ^ K+iy9) < 1.3 -10-^ [55], (123) 

both at 90% C.L. and assuming the SM differential rates. 

However, in general, the kinematical distributions (and the associated phase-space ranges) depend on the 
nature and couplings of the invisible particles. So, one would need to correct for the associated experimental 
reconstruction efficiencies and background shapes. Most notably, the SM backgrounds typically rise steeply at 
small final state K^*^ momentum in the center of mass frame, thus reducing signal sensitivity in this region [82] . 
Without the detailed knowledge of the experimental analyses and detectors, we cannot faithfully reproduce the 
final signal sensitivity distribution. However, we consider these effects to be the least severe for massless final 
state invisible particles, since the kinematical distributions in this case at least cover the whole kinematical 
region populated by the SM signal with (almost) massless neutrinos. There we derive our tentative bounds on 
the individual NP operators given in Tables |4| and |6| The impact of purely kinematical (phase-space) effects 
on the NP scale sensitivity away from the massless invisible particle limit is illustrated in Figure [2| for the case 
of pair production of two invisible fermions through the various considered operators. 

The Super-B factories are expected to provide a sensitivity down to a fraction of the expected SM signal. 
Then, possible measurement of the K^*"^ momentum or missing energy distributions could provide additional 
powerful discriminants in the search for NP contributions f5^. In this respect, the decay B — > K*X has 
the virtue that the angular distribution of the K* decay products allows to extract information about the 
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polarization of the K* . The experimental information that can be obtained from the process B — K*{-^ Ktt)X 
with an on-shell K* is completely described by the double differential decay distribution in terms of the two 
kinematical variables s = {pb — Pk)^, corresponding to the invariant mass of the final state invisible particles, 
and 9, the angle between the K* flight direction in the B rest frame and the K flight direction in the Ktt (K*) 
rest frame. The spectrum can be expressed in terms oi B K* transversity rates T^ t corresponding to 
longitudinally and transversely polarized K* final states (see App. lB.l~4l) . while the double differential spectrum 
can be written as 

— ^ -^siTiH+-—^cos^9. 124 

ds d cos 9 4 ds 2 ds 

Thus, dTL/ds and dTx/ds can be extracted by an angular analysis of the K* decay products. Finally, the total 
invisible mass distribution is 

dV , „ d^T dVr dTL 



, / dcos9— n^-r + -r^- (125) 

ds J_i dsd cos ds ds 

As an illustration of the potential impact of such future precision measurements, in Tables E] and El we also 
present the accessible NP scales assuming a 20% relative precision on the B — )• K^*^X rates compared to their 
SM predictions. 



B.1.2 b ^ d+ missing energy 

The B^ — > 7r+(p+)X processes in the SM, with X = i>D, receive a dominant contribution already at the tree 
level though the decay chain _B+ r^Vr 'K^{p'^)vrVT [HI]- They have recently been measured by both 
Belle IT and Babar [16]. Reinterpreting their B{B^ — t^Vt) values in the 7r(p) decay channels of the r by 
multiplying them with the corresponding well measured r — >■ tt{p)v branching fractions we obtain 

B{B+ ^ Ti+vv) = 1.96(85) • 10"^ , (126) 
B{B+ p+vv) = 0.97(50) • lO""* . (127) 

A potential NP signal in these modes would manifest itself via differences in the measured B{B ^ tv) 
values in the leptonic and hadronic decay modes of the r - something that future dedicated experimental 
searches could employ to reduce systematic uncertainties. Although the SM signal shape in this case is well 
determined and largely free from theoretical form factor uncertainties, the appearance of two neutrinos in 
the final state means that the same experimental caveats in extracting bounds on NP contributions apply as 
in the B — > K''*^X case [84|. Again the polarization states of the p in the B — > pX mode, which are well 
predicted within the SM, could be reconstructed using the angular distributions of the p — >■ 27r system and aid 
in discriminating SM contributions from possible NP effects. 

The neutral modes B^ — 7r°(p°)X are free from long distance SM contributions, but the purely neutral 
final states make them more challenging experimentally. The present bounds of 

B(S° ^ 7rVj>) < 2.2 • 10"^ [H3] , (128) 
B(B° ^ p^vD) < 4.4 • 10^^ [83] , (129) 

are less constraining than the charged modes analysis. Consequently, in setting our bounds on invisible particles 
in the massless limit in Tables Eland [5] and away from this limit in Figure[51we tentatively allow NP contributions 
to saturate the experimental uncertainties in the charged modes, Eq. (I126p . 



B.1.3 Other modes with missing energy 

Both Belle [SI] and Babar [5B^ have searched for the B — > XX decay mode, which is helicity suppressed in 
the SM. While unresolved soft photons can partially lift this suppression [HIIHZ]: the SM predictions for the 
branching ratio remain at the order of 10~^ ;6j. Being a two body decay process in the scenarios we consider. 
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the kinematical distributions are trivial and no model-dependent efficiency corrections are needed. The latest 
experimental upper limit reads 



B{B ^ XX) < 1.3 -10-^ @90% C.L. |85] , (130) 

and can be employed directly to constrain the relevant interactions of invisible particles, as given in Tables S] 
and [H and Fig. [51 In the future, both Bg^d — > XX modes could potentially be probed to greater accuracy at 
the super flavor factories [88] . 

On the other hand, while Bs^d — > jX with reconstructed final state photons allow to lift the helicity 
suppression suffered by the NP fermionic contributions proportional to Jaa, the additional aEM suppression 
only makes them competitive with the B — p{K*)X modes in the opposite region of large invisible particle 
masses (m^ — mp(^x-') < nT-x < itt-b)-, where the helicity suppression is least and the Bs^d — ^ X channels are 
effective in constraining NP effects as well. In addition, precisely in this region the SM backgrounds, most 
notably from misreconstructed Bq — >■ ^X events, limit the experimental reach of i?s,d "fX [86| . For a more 
detailed discussion of the expected NP sensitivity, we refer to Ref. [6 and do not consider these modes any 
further. 

B.1.4 Transversity basis in B H* transitions 

Consider the decay B — )• H*X with the B meson decaying to an on-shell vector meson H* and other (invisible) 
particles X. The amplitude for this process can be written as 

7W(„)(B^i/*X) = e^4m)Af^, (131) 

where e^. (m) is the polarization vector of the H* . Being on-shell, the H* has only three polarization states, 
satisfying e_f/. ■ ph* = 0. In addition, the polarization vectors satisfy the following relations 

e*j^. {m)eH'tj.im') = -Smm' , ^ e^. i^^)<^H' i'm')6,nm' = -g^" + ^^'f^* • (132) 

m,m' ^ 

In measurements where only the i/* decay products are reconstructed, only two of the three i/* polarization 
states can be disentangled. The corresponding transversity rates {Tl and Ft) can be projected using covariant 
polarization projectors. First, we write the sum of squared amplitudes over polarizations as 

\M(^ra){B H*X)\^ = M^MlY^"' , P^'' = -y^'' + ^^"f^" . (133) 

m -H 

The longitudinal H* polarization vector should satisfy p//* \ \tH* in the B meson rest frame and the correspond- 
ing decay rate (F^) is obtained with the help of the projector 

( ^ PB ■ PH* ^ \ ( PB ■ PH- 



where niBy^Jj^^l'^ — |Pif|, \z = A(l,2;,rf), is the absolute value of the B* momentum in the B rest frame. 
The (unpolarized) transverse B* helicity projector entering the transverse rate Ft is then simply obtained as 

B.2 Matrix elements for B decays 

The hadronic matrix elements between a Bq state and the vacuum are 

(0|67M759|-Bq(pB)) = , (135a) 

(0|575g|B,(Pi3)) = f^/B, . (135b) 

(mb -I- mq) 
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For the decay constants, we use the values of a recent lattice QCD average [SI]: fs^ = 0.2388(95) GeV and 
fsa = 0.1928(99) GeV. The matrix elements for B ^ Hq (where iJ is a pseudoscalar meson, i.e. K or tt) 
are [SO] 



{H{pH)\b^y\B{pB)) = n{r) 



+ /o^(9^)^^9^ 



{H{pH)\bq\B{pB)) - 

{H{pH)\ba^'''q\B{pB))=^ 



rrib - ruq 



m_B(l + th] 



(136a) 
(136b) 
(136c) 



where — {pb+phY ^^'^ 9^' = {pb^PhY ^ while z — q-^ /m\ and r; — mi/niB- The {H [pH)\b(T^^^ ^^q\B{pB)) 
matrix element can be obtained via the Chisholm identity. We also consider B H* matrix elements [OT] 
(where H* is a vector meson, i.e. K* or p) 



{H*{pH^,eH')\bj^q\BipB)) = 2e^,p,ejr*P>^ 



"^b(1 + r/f.) 



(i/*(p^^.,eff.)|fe7'^75g|i?(ps)) =i4''*"iB(l + rK)Af {q') - zP'^ie*^. ■ q) 



mB{l + Tif. 



^q^^{e*B,■q)^-^[Af{q')-Af{q')], 
rnsz 



{H*{pH.,eH')\h^q\B{pB)) ^ -2i- 



niB 



rK{e*H. ■ q)A^' [q^) 



nib + niq 

{H*{pH*,eH*)\ba^^q\B{pB)) = ie^^^pa 



pPe*if. 



z 



qPP" 



Tfiq') 



. 1 — r 



(137a) 

(137b) 
(137c) 

(137d) 



while the matrix elements {H*{pH* ,eH*)\ba'^'^q\B{pB)) are again determined via the Chisholm identity. At 
q'^ = the form factors satisfy the following relations /f (0) = /(f (0), ^^'(0) = (0), and T/^* (0) = 
(0) — (0)- For convenience we will also define the following auxiliary form factor combinations 



iq')^A^ iq-)-^il + rH.ril-ri,.~z)A^ (g^) , 



1 



Tf (q^) ^ ff [q') ^(1 - rl.){l - rj,, - z)Tf{q^) 



(138a) 
(138b) 



where Xiz = A(l, z, rf). 

In our numerical analysis, we employ the form-factor normalizations, shapes, and the associated uncertain- 
ties as determined in Refs. [501 - 1^ using light-cone QCD sum rules. In particular we employ results of Ref. [HI] 
for the V^''P,Af 'P,Ti values, Ref. [90] for the /f^o ^ and form-factors, while we use the results of a 
more recent calculation [32] for q. 
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B.3 Spin 1/2 invisible particles in the final states 

The differential rates into a pair of invisible fermions are 



+j;^'''^mAvf^r) + Jii'^^'^'^'mspf^A) + jr^^i/^^i' + . (i39d) 

where vf" = /|/i/87rmB„ rf = mBl3i,X]^f,,j96n^ A. = ' "^rf/z, A = y/T^^ and 

= <Ji^b + m,f ,1^^ = X^^'pl ,1^^'^^ = 4r^mB,/(m6 + m,) ,1^^ = 4^ , 
4 = f^{q'fzXH.l3lJ2{l + r„f , = j/(l + Sr^, /z)/^^^ , J^""'^ = 3/f (92)/+(92)r-^AH./(l + vh) , 
JV = Ai/z/f (9^)^(1 + 2r|/.)/2 , = jyfi%J{\ + 2ry.) + 3/o^(g^) V^(l - rlfjz , 
J^"- = 3^/0^(9^)2(1 - rlfml/A{m, - m,f , = J^^^^^ , j;f ^'^^ = J^''Ar^{m, - m,)/mBZ , 

= Al,.,Tf *(g2)2(i + Srl/z)/Ar%.z,j:p^ = X]j,Jl;Tf{q^f/Ar%.z, 
J'f''^ = 3r^X%,,Af {q^)Tf {q^)/2r%, (1 + th^ )z , ^ = SzXn.Af {q^fm%/4{mt + m,f , 
JT = ^r^l^ ' ^r'^^ = ^H^H-.Af {q^fmB/2{mB + m,) , 

^^'^ = Al„,Af {zf{l + 2tIIz)I%t\, (1 + TH' f , j;^^ = 4^^^^./(l + 2ry^) + 3AH..r2 {q'f/z , 
Jf = 2XH^^PlTr{qr + (1 - r^^rrnqn^ + 8rV^) , ^r'' = AH..F^*(g^)(. + 2r2)/(l + r^O' , 
4'^ = 2AH..Tr (9^)^(1 + 8rl/z) + (1 - r|,0'/3i7^2^*(9^)^^^^'^ = 12r^AH..Tf *(,2)y(,2)/(i + ,^.) , 

= AH.../3i^^*(g^)/(l + rH.)^^^^■^ = 12(1 - r%,){l + r„.)r^Ar{q')Tr{q') , 
j;^^ = (1 + rn. )'^Ar (g2)2(i + 2^/^) , J^'^^ = J^'^^^^ J(l + 2rl/z) . (140) 



B.3.1 Standard Model rates 

The -B+ — > 7r"'"(p+)i^i/ modes are dominated by the tree- level contributions mediated by an intermediate on-shell 
tau lepton. To an excellent (small tau-width) approximation, they arc given by 



^{B ^ H*i:{p) 




(139b) 



(139c) 




— (B+ w+ui?) 



dz 



(141b) 



(141a) 



(141c) 
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where Fgj^'''' = "^B?'r/s/^(p)l^''Vl^(^pG^/647^^^^-. To obtain the total rates one needs to integrate over the 
available phase-space (z), which for the on-shell tau contributions is given by z G [0, (1— r'^)(l — r^j-p^/r^)]. Alter- 
natively one can normalize these distributions to the experimentally determined B ^ tv (from purely leptonic 
tau reconstruction) and r — T^{p)v (from prompt tau decays) branching ratios since B{B^ — !■ Tr^{p^)v9) ~ 
B{B — > Tv) X B{t — > T:{p)iy), eliminating the theoretical uncertainties in the normalization factor Fq. The 
neutral modes are dominated by short distance loop contributions, similar to i? K^*^i>D, albeit further 
CKM suppressed, leading to branching ratios almost two orders of magnitude smaller compared to the charged 
modes [84j . 

The dominant (short distance loop) contributions to the kaon modes {B K^*^i>D) are given by 

^{B ^ K.D) = Ti^xT.Uiq'r , (142a) 



'^^T/r^ r^i, _N T^BKn ^ 1/2 

— {B ^ K*vv) = V^^2z\j^.^ 



(l + r^.)' 



(142b) 



where now T^^ = m%{GFa\VtbVt%\Clff/2bQT:'^ and \Clf\ = 6.33 ± 0.06 [7i[93]. The charged modes also 
receive tree-level contributions mediated by intermediate on-shell r leptons similar to i?+ — > tt'^{p'^)vv. These 
are however always subleading in the kaon case; of about 15% of the above short-distance contributions |84) . 

B.4 Spin invisible particles in the final states 

The rates for the production of a single invisible scalar from the H* {DQ^ (p operator of Eq. (PD|) are 

T{B ^ H<t>) = ^i"\9s? l^^~"'fy ^Hf , (143a) 
Tl{B ^ H*<t>) = rf/'igpp ,, Af{^lf , (143b) 

where Ffj-'^**' = msA^;^,) ./167r, \ij^z = A(l,rf ,r^). The vector operator contributions are related to these via 
quark EOM, see Eq. The rates for the production of two scalars are 



'Tip „ kpsPn^R 

^iB~,H^<P)^^T^,"(^J,^%ss\' + Jr^\9vv\'} . (144b) 
^{B^ H*U) = ^T^' [jr\9PS? + Jr^\9AV? - Jr''^^{9Av9*Ps)] , (144c) 

where 

J^' = J^^iz^jJ^ = \H.PlJ^{q'fl^,Jr = Jri^^^^r = ^r^/30./4(l + ir^^lz) . (145) 
Note that if (/> = i^, Bose statistics has to be enforced and these rates should be divided by two. 
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B.5 Spin 1 invisible particles in the final states 

The production of a single vector using the simple FCNC operators of Eq. ([25)) or the gauge-invariant operators 
of Eq. dMl) are 

B^HV: rPl = rfy"jrv\ev\' , r™ = !I^Tfy"jT^\hT\' , 

B^H*V: r^^=Tfy^'j;^\eA\\ rf^ = ^r?v"'Jit\hf?, (146) 

2 

_ pSH* r -T«y I, |2 I -T"A|, .|2\ _ ^B tBH' f ^r/T\u |2 I 'T//T| L I2\ 

^T-'^iv Wivl'^vl +Jiv\(-a\ ) , It ""^xT^iv +Jiv\'T't\ / ' 



where 



^/'^ = 2(1 + r^O'^r("^y)' = 44Aw/#(™2y) V(l + ^v)' , (147b) 

Jlf = Al^.^Tf* (m2^)Vr^r^. , Jl'^ = SA^.v^Tf* (rr4)^ = 8(1 - rj^.^Tr {mlf . (147c) 

Note that in the my — >■ limit, F^^' (i? ^> H*V) is well-defined while the tensor operator contributions actually 
vanish thanks to the form factor relation T^(0) = 0. 

The rate and differential rates for the production of two vectors are 

r["l(i?, ^ VV) = ^r^^ {ir\hpp\' +I^'\hps\'} , (148a) 

^{B ^ HVV) = {Jr\hsp? + J^'\hss?] , (148b) 



"^'^ {B^H*VV)='^Tr' {J{r'\hps? + J{r''\hpp\^] . (148c) 



where 



dz ^ ' 



tIPS o2 ^IPP ^IPP rriPP 



jr = zp'v.jr ' = ^^jr^Pv. + 64/^^) • (149) 



B.6 Spin 3/2 invisible particles in the final states 

The rate for the '5^' modes are 



T{B, ^ m ^ ^rj' {l7|/pp|2 +I^,^|/P5p -i;,-^-^^5i(/^^/;p) , (150a) 

^(5 ^ i/*^) = !^rr {4\ff^? + Jl\hTf + Jj,'\fTs\' + 4,'\fts\' 

rTPif |2 I ^fP\r^ |2 I ^VV\r |2 , ^VA\ r |2 , ^SP\r |2 I -7SS1 



I-^^-pI I'^tpI +'7'^' +J'^/ \fvA\ +J^,' \fsp\ +J4,' \fss\ 



-^.''-''''mfpfvv) + jV^mvviTT) + j^^^'^mvArsp)} , (isob) 
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where I3[ = (9 - Q^j + 5/34)/18, /3f = (5 - 6/3| + Qpt)/lS and 

j} = z^jjl3'^j4 ,1 = z^J^il + 2rl/z - Ur^/z-" + 80r|/^3)/9/3|, , = + 5ry^)^2^/36 , 

= zj^rlii + lOr^ /z)/36 . J^^ = z^J^{l + 3r|/^)/36 , jj^ = 5^ jjr|(l + 2rl/z)/^&0l, , 
jy = Xn.z-'fl {q^){mrl/z^ - 2r%lz^ - 2rl/z + 1)/18 , 

jy = z[M{tI - \fr%!^{qffi'^^ + \H.zfl{q^ffi%^{^fi'^^ - z)]/12 , J^F = z^J^P'^/^, 

J'/ = ziJ'^Pi^;j4 ,ir^^' = 2zy^/9 , jp^' = zy^i - ior|/.)/9 , j^' = ^i? , jir = ^jy , 

= z'jjil + rl/z)/9 , Jf'^^ = jr^^z^P%j21 , 4^'^^ = Jj'^'^z^l + 2rl/z + 6r^/.2)/54 , 
J^'^ = Jr^z\Z - Url/z + 38rV.^)/27 , J^'^ = jf'''^z%jA . (151) 

With the many possible interferences among the tensor currents and the comphcated kinematical functions, 
the differential rates for B H*^^ are too cumbersome to be given here. 
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